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Hadronic form factors for the rare weak transitions Af, — >■ A'*' are calculated using a nonrela- 
tivistic quark model. The form factors are extracted in two ways. An analytic extraction using 
single component wave functions (SCA) with the quark current being reduced to its nonrelativistic 
Pauli form is employed in the first method. In the second method, the form factors are extracted 
numerically using the full quark model wave function (MCN) with the full relativistic form of the 
quark current. Although there are differences between the two sets of form factors, both sets satisfy 
the relationships expected from the heavy quark effective theory (HQET). Differential decay rates, 
branching ratios and forward-backward asymmetries (FBAs) are calculated for the dileptonic decays 
Aj, — >■ A''*'^+£~, for transitions to both ground state and excited daughter baryons. Inclusion of 
the long distance contributions from charmonium resonances significantly enhances the decay rates. 
In the MCN model the A(1600) mode is the dominant mode in the /i channel when charmonium 
resonances are considered; the A(1520) mode is also found to have a comparable branching ratio to 
that of the ground state in the ^ channel. 

I. INTRODUCTION 

The weak decays of heavy hadrons have been an important source of information on some of the fundamental 
parameters of the Standard Model (SM). In particular, semileptonic decays of heavy hadrons have been important 
in the extraction of some Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [1]. However, the description of 
these processes requires a number of a priori unknown form factors that parametrize the uncalculable (to date) 
nonperturbative QCD dynamics. The precision to which these form factors can be calculated or modeled limits the 
accuracy with which the CKM matrix elements can be extracted. In this regard, the heavy quark effective theory 
(HQET) [2], quark models [3, 4], QCD sum rules [5], lattice QCD [6], etc. have been employed to improve the 
modeling of these form factors. 

In the same way, decays of heavy hadrons induced by flavor-changing neutral currents (FCNC), the so-called rare 
decays, have been a subject of significant interest in recent years. Processes like b ^ sj and b — s£'^i~ are forbidden 
at tree level, and at one-loop level are suppressed by the Glashow-Iliopoulos-Maiani (GIM) mechanism. These decays 
therefore receive their main contribution from one-loop diagrams with a virtual top quark and a W boson. Therefore, 
they provide valuable information about the CKM matrix elements Vts and Vtt- 

In addition, because these decays occur at loop level, they are sensitive to new physics beyond the SM. In these rare 
decays, new physics can appear either through new contributions to the Wilson coefficients that enter into the effective 
Hamiltonian that describes these decays, or through new operators in the effective Hamiltonian that arise from sources 
beyond the SM. For these reasons, these decays are promising candidates for looking for new physics beyond the SM. 
Thus, the values of the Wilson coefficients are crucial to the determination of any new physics. As with the case of 
semileptonic decays, the accuracy to which these parameters can be extracted is limited by our knowledge of the form 
factors that are used to parametrize the hadronic matrix elements. Many experimentally measurable quantities such 
as branching ratios, forward-backward asymmetries, lepton polarization asymmetries, etc. can be analyzed to more 
precisely determine SM parameters and to look for new physics beyond the SM. It is known that most observables will 
depend sensitively on the form factors; thus, the accuracy to which these form factors can be calculated is paramount 
to the determination of any new physics inferred from these decays. 

There have been many theoretical investigations of rare 5 — s transitions in the the meson sector [7-52] . Melikhov 
et al. have used a relativistic constituent quark model to treat the dileptonic decays of the B meson [21]. In [41], 
perturbative QCD (pQCD) is used to estimate the form factors for B — K'^*K Light-cone sum rules (LCSRs) have 
been employed for these transitions as well [15, 29, 35, 46]. The form factors obtained from the various models 
have been used to calculate observables such as branching ratios (BRs), forward-backward asymmetries (FBAs) and 
lepton polarization asymmetries (LPAs). These observables have been calculated both within the SM and for various 
scenarios that arise beyond the standard model, such as supersymmetric (SUSY) models [29, 34, 42], models with 
universal extra dimensions (UEDs) [52], and other new physics (NP) scenarios [44, 46]. 

In the baryon sector, there have been a few techniques employed in the computation of the form factors for the rare 
Af, A transitions [53-62] . In much of what has been put forth, HQET is used to reduce the number of independent 
form factors to two universal form factors valid for all currents, and a model is then employed to compute these two 
form factors. Aslam et al. [53], Wang et al. [54] and Huang et al. [55] have used LCSRs to obtain the form factors for 
Af, — > A. QCD sum rules (QCDSRs) have also been employed in obtaining these form factors [56-59, 62]. A number 
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of authors [56, 58, 60] have also used pole model parametrizations (PM) for the form factors. In [61] the form factors 
for these transitions have been estimated using pQCD. The MIT bag model (BM) has been used in [59]. To the best 
of our knowledge, no transitions to excited state A's have been explored. 

Experimentally, the exclusive radiative transition B — >• if* (892)7 was first observed by CLEO [63]. The branching 
ratio for this mode has been measured recently by both Belle [64] and BaBar [65] with an average branching ratio 
B{B° is:0*(892)7) = (4.33 ± 1.9) x IQ-^. Several other decay modes such as B ^ Kij, K^{U30)j, etc. [66-69] 
have been observed as well. The mode JB° — > <j)j has been observed by Belle; its branching ratio was measured to 
be B{Bg — > (^7) = (57^19) X 10~^ [70]. Branching ratios for the inclusive process B Xg^ have been measured by 
CLEO [71] , BaBar [72] , and Belle [73] . Prom these measurements, the average branching ratio is [74] 

BiB -s- X,-f) = (3.55 ± 0.24 ± 0.09) x 10"*. 

Branching ratios for the dileptonic decays have been measured by Belle [75], BaBar [76], and CDF [77]. Prom the 
BaBar and Belle data, B{B Kl+i') = (0.45 ± 0.04) x lO'^ and B{B K*{892)e+£-) = (1.08 ± 0.11) x 10"^ 
Branching ratios for the inclusive process B — > Xsi~^i~ have also been measured by both Belle [78] and BaBar [79]; 
the average branching ratio is [80] 

BiB ^ XJ+e-) = (3.66t°;^^) x lO"*^. 

Note that this means that decays to the K and K* account for less than 50% of the rare dileptonic decays of the 
B meson. The Belle Collaboration [81] has also measured the forward-backward asymmetry in B — )• K*i~^£~ and 
extracted ratios of Wilson coefficients from those data. 

The experimental situation for 6 — ^ s transitions in the baryon sector is less rich than in the meson sector. The CDF 
collaboration recently reported the first observation of the baryonic PCNC process A;, A/i+/U~ [82]. The branching 
ratio for this mode was measured to be B{Ab A/x+^") = (1.73 ± 0.42 ± 0.55) x 10~^. The LHCb Collaboration 
[83] estimate that with 2 fb~^ of data taken in one year, there should be 750 A(, — >■ A7 events and nearly six times as 
many events in excited hyperons. For the dileptonic decay mode, 800 events are expected for decays to the ground 
state A; no potential yields were reported for excited states. 

In this paper, we examine the rare weak dileptonic decays of Afc baryons to ground state and excited A baryons. 
Some of the motivation for this study has been outlined above. It will also be useful to examine the sensitivity 
of some of the lepton asymmetries to the form factors. It has been shown for rare meson decays that there are 
observables, such as lepton polarization asymmetries, that are largely independent of the form factors in certain limits 
[7, 12, 13, 15 18]. Such quantities thus offer largely model-independent ways to examine the physics content of some 
of the Wilson coefficients [13]. 

To this end, we use two approximations to compute the form factors for Aft — )• A^*' transitions. First, we use the 
approximation employed in [84, 85]. In that work, analytic form factors for Ab Ai*\ Ac A^*^ Ah A^(*) and 
Ac — >■ A^(*) were calculated using single component wave functions obtained from a variational diagonalization of 
a quark model Hamiltonian. In that calculation, quark operators were reduced to their nonrelativistic Pauli form. 
The decay rates obtained using the form factors extracted in that calculation were in reasonable agreement with 
experimental results for the semileptonic decays Ah AciUi and Ac Ktut . 

The second approximation we use for computing the form factors is an extension of this method; we keep the full 
relativistic form of the quark spinors and use the full quark model wave function in a numerical extraction of the 
form factors. Because this method uses fewer approximations, the form factors obtained in this way should give more 
reliable results than the form factors obtained from the first method. Wc present the results of both methods to 
demonstrate the sensitivity of the observables to the form factors. The goal here is to find observables that may be 
less dependent on the form factors. 

The rest of this paper is organized as follows: in Section II, we discuss the hadronic matrix elements, decay rates, 
and forward-backward asymmetries. Section III gives a brief overview of HQET and presents HQET predictions for 
the relationships among the form factors for the transitions we investigate. In Section IV, we describe the quark 
model used to obtain the form factors, including some description of the Hamiltonian. The two methods we employ 
to obtain the form factors are also briefly discussed. Numerical results such as form factors, differential decay rates, 
branching ratios, and forward-backward asymmetries are presented in Section V. We present our conclusions and 
outlook in Section VI. Some details of the calculation are shown in the Appendices. 
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II. MATRIX ELEMENTS, DECAY RATES, AND FORWARD-BACKWARD ASYMMETRIES 

A. Matrix Elements 

The amplitude for the dileptonic decay of the A;, baryon can be written 

iM{A, ^ Ain-) = ^^VMH^L<y^ + H^Llf^), (1) 

where iJf and H2 contain the hadronic matrix elements 

= -^C7(m6)Tj^ + C9(m6)J^, (2) 
= (7io(m(,)J^ (3) 

The Ci are the Wilson coefficients, 

= (A(pa,sa) I +75)6 | Ab(pA,,SAj>, (4) 

and Jl is the matrix element of the standard V — A current 

= (A(PA, Sa) I S7'*(l - 75)6 I AbipA^SA,)). (5) 

L^^^ and L^jJ^^ are the vector and axial vector Icptonic currents, respectively, written as 

L*,^) = ue{p-,s^hf,ve{P+,s+), (6) 
^If^ = ue{p-,s~)-fij,j5ve{p+,s+), (7) 

where arc the spin projections for the Icpton and the antilcpton. 

In our analysis of the dileptonic decays, we will include the long distance contributions coming from the charmonium 
resonances J/ip, ip', ... etc. To include these resonant contributions, we replace the Wilson coefficient Cg in Eq. 2 
with the effective coefficient 

= Cg+ Ysd{z, s') + Yld{s'), (8) 

where z = mc/mi, and s' = q^/ml. YgD contains the short distance (SD) contributions from the four-quark operators 

far from the charmonium resonance regions and y^/j are the long distance (LD) contributions from the four-quark 
operators near the resonances. The SD term can be calculated reliably in the perturbative theory, but the same cannot 
be done for the LD contributions. The LD contributions are usually parametrized using a Breit-Wigner formalism by 
making use of vector meson dominance (VMD) and the factorization approximation (FA). The explicit expressions 
for YsD and Yld are [7, 53, 54, 56-58] 

Ysd{z,s') = (3Ci+C2 + 3C3 + C4 + 3C5 + C6)/i(z,.s')- 

i(4C3 + 4C4 + 3C5 + C^)h{l, s') - ^(Cs + 3C4)/i(0, s') + 

^(3C3 + C4 + 3C75 + C6), (9) 
Yld{s') = ^(3Ci+C2+3(73 + C4+3C5 + C6) 

(10) 
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where mj, Fj, and r{j — )• £~^£ ) are the masses, total widths, and partial widths of the resonances, respectively, 
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and kj are phenomenological parameters introduced to compensate for VMD and FA; they are chosen so as to 
reproduce the correct branching ratio of B{B K^*^V K^*^f+e-) = B(B K'-*^V)B(V f+f~), for V = 
J/ip, 4''- Since none of the analogous branching ratios of the A;, have yet been measured, we apply the phenomenological 
factors obtained from decays of the B mesons to the decays of the A;,. For the lowest resonances J/ip and tp', we use 
k = 1.65 and k = 2.36, respectively; for the higher resonances, we use the average of J/tp and tj/ (see [35, 62]). 

The matrix elements in Eqs. 4, and 5 contain the currents sj'^b, 57^756, sia^'^b and 'sia^'^^^h. In this work, we 
examine decays to daughter baryons with = 1/2+, l/2~, 3/2~, 3/2+, 5/2+. We make this choice because in the 
quark model we use to calculate the matrix elements, these states have the most significant overlaps with the initial 
ground state within the spectator quark approximation. 

A baryon with angular momentum and parity may be represented by a generalized Rarita-Schwinger field (or 
spinor-tensor) Uy^i.../j„ (pa), with n = J — 1/2 indices. These spinor-tensors are symmetric in all of the indices, and 
satisfy the conditions 



|^a'«,.i...m„(pa) = toa'«^,...^„(pa), 7''''"^i.../x„(Pa) = 0, 



Pa «w-m»(Pa) = 0. 5'' 



*'^l/...(X- 



.(pa)=0. 



(12) 



For parity considerations, it is necessary to divide the spinor-tensors into two classes. Those with natural parity are 
called tensor, while those with unnatural parity are labeled as pseudo-tensor. Here, a state of total angular momentum 
J is said to have natural parity if P = (— l)'^"^/^^ unnatural parity otherwise. 
For transitions between the ground state and any state with = 1/2+ , the matrix elements of these currents are 



(A I s7''6 I Kb) = u{pA, sa) 
(A I S7''756 I Afe) = u(pA, sa) 



Fi(g2)7/' + ^2(9')^;'* + ^^3(9')^^'" 



«(PAi,,SAj, 

75u(PAb,SAj, 



(A I 52(7^-^6 I Ab) = u(pA,SA)r^"M(pA,,SAj, 

where we have used v = Phb/'mAti ''j' = Pa/'t^Aj and 



Hi{q^){v''v"' -v^v'*"). 



(13) 



(14) 
(15) 



(16) 



The Fi, Gi, and Hi are the form factors of interest, and these are functions of the square of the four-momentum 
transfer q^ = (p^^ — pa)^ between the initial and final baryons. Since 



(17) 



the matrix elements involving the current sia'^'^^sb can be related to those involving sia^'^b. 
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Similarly, the matrix elements for decays to daughter baryons with = 3/2 are given by 



(A I 57^6 I Afe) = u„(pa,sa) 



Fa9 



OtjJ, 



m(pa,,saJ, 



(A I S7''756 I Afc) = u„{pk,sa) 



75M(PA6,SaJ, 



(A I sia'^'^fe I Ab) = m„(pa,sa)T"^"«(pa„saJ, 



where 



For decays to = 5/2+, the matrix elements are 

(A|S7^6|A6) = Ua0{PA,SA)v°' 



F17'' + F2V^' + F^v'^ 



Fig^" 



•"(PAmSaJ, 



(A I 57^756 I Afe) = Ma/3 (pa, Sa)!"" 



yf" (^Gi7'^ + G2f ^ + Gaf '^^ + 



G4/'' 



75'w(pa,,saJ, 



(A I SZa^^H A6) = Wa/3(PA,SA)r"^''''w(pA„SAj, 



where 



(18) 



(19) 
(20) 



(21) 



(22) 



(23) 
(24) 



(25) 



Thus far, only the matrix elements involving spinors with natural parity, i.e. spinors with parity (— 1)'^^^^^, have 
been presented. The equations involving spinors with unnatural parity can be found by inserting 75 to the left of the 
parent baryon spinor in the equations for natural parity. 

The matrix elements for the tensor currents can be written in a more convenient form. If Eq. 15 is contracted on 
both sides with the four-momentum transfer q^, use of the equations of motion leads to 



(A I sia'^'^q^b \ Ab) = u{pA, sa) 



Fl{q^)r + F^{q^)v^ + Fi{q^y 



w(pa6,saJ, 



(26) 



where, 



= - (n^Ab + "^a) Hi - {niA, - niAV ■ v')H2 - (toa,," • v' - mA)Hs, 
F2 = ruA^Hi + (mA, - toa) H2 + {ttia^v ■ v' - mA)i?4, 
Fg^ = tuaHi + (rriAfc - ttia) -ffa - (wa^ - niAV ■ v')Hi, 



(27) 
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valid for states with = 1/2+ . For the axial tensor, a similar procedure leads to 



(A I sia'^^j^q.^b \ Aj,} = u{pa, sa) 



75M(PAfc,SAj, 



(28) 



with 



Gj = {rriA,, - ruA) Hi - mA{l - V ■ v')H2 - mAi,{l - V ■ v')H3, 
G2 = ruA^Hi - mAH2 - mA^Fg, 
G3 = rriA-ffi + mAi?2 + mA„Hs. 



(29) 



Similarly, for transitions to states with = 1/2 , the matrix elements for the tensor and axial tensor currents 
become 



(A I sia'^"q„b \ A^) = u{pa,sa) 
{A\sia^'''y5qt,b\Ab) = u{pa,sa) 



Fl{q')r + Fl{q')v'^ + Fl{q')v'^ 

Gl{q'')r + Gnq')v'^ + G^iq^)v'^ 



75W(pAi,,SAj, 

U{PA,,SA,), 



respectively, where 



and 



Fi = (mAi. - toa) Hi - (ttia, - ttiav ■ v')H2 - (niA^v ■ v' - mA)Hs, 
F2 = mA^Hi - (mAf, + m.A) i?2 + {mA^v ■ v' - mA)Hi, 
F[ = uiaHi - (mAt + itia) H3 - (rriA^ - uiav ■ v')Hi, 



G1 = ~ (mAt + mA) Hi + mAi^ + v ■ v')H2 + mAh{l + V ■ v')H3, 
G2 = mA^Hi - TOAi?2 - niAkHs, 



G3 ttiaHi ~ mAi?2 - WAt-ffa- 
For transitions to states with = 3/2~, we obtain 

{A{pa) \ sia'^" qt,b \ Ab{pAi)) = «a(PA,SA) 



Fjg'^^ 



w(pa,,saJ, 



(A(pa) i sia'^"^5q^b \ AbipA^)) = UaipA,SA) 



G'^j^ + G'^v^ + Gjv'^ ] + 



75'"(PA6,SaJ, 



while for transitions to states with = 5/2+, we have 

(A(pa) I sia^^q^b \ AbipA^)) = Ua0{PA,SA)v°' 



■«(pa6,saJ, 



(A(pa) I sia'^'^^sq^b \ AbipA^)) = Ua0{pA,SA)i 



v^(gJ^'' + GJv^ + GJv'^^ + 



75M(PAfc,SAj, 



(30) 
(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 
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where Ff and Gj are given by 



FI = - {niA^ + rriA) Hi - {rriA^ - m^v ■ v')H2 - {niA^v ■ v' - mf^)Hs - rriA^H^, 

F2 = mf,,Hi + (rriAi, - tua) H2 + {mA^v ■ v' - mA)Hi - niA^He, 

F3 = mAHi + (mA, - ttia) H3 - {niAt - ^av ■ v')Hi, 

FJ = {mAi - mA) H5 + (niAt, - mAV ■ v')He, 

Gf = (mAfc - ttia) Hi - mA(l - v ■ v')H2 - tua^ {1-v- v')H3 + mA^-ffs + mA-ffe, 
Gl 

Gl 
Gl 



mAfc-ffi - mA-ff2 - mAfc-ffs, 

rriA-ffi + TOA-ff2 + m\, Hj, - rriA-ffe, 
(toa6 + toa) -ffs + wa(1 + V • v')H^. 



For = 3/2+, a state with unnatural parity, we have 

(A(pa) I 'sia^^qj) \ Ab{pAj) = Ua(PA,SA) 



v'^lF^jf^ + F^v'^+F^v"' 



(g^j" + G^yf" + Gjv''' 



Glg'^i' 



u(pAfc,SAj, 



with 



FI = (mA(, - toa) Hi - {mA^ - mAV ■ v')H2 - {m^A^v ■ v' - mA)H3 - mA^H^-, 

F2 = mAtHi - (mAfc + niA) H2 + {niAtV ■ v' - mA)i^4 - WA^-ffe, 

F3 = TUaHi - (piAi + TflA) Hs - (jTlAfc - niAV ■ v')H4, 

FJ = - (rriAi, + ruA) H5 + {ruA^ - niAV ■ v')Hq, 

GJ = - (mA(. + ttia) Hi + mA(l + v ■ v')H2 + mA(,(l + v ■ v')H3 + mA,,H^ + mA-^e, 

G2 = mA^Hi - mAi?2 - mA.Hs, 

GT = vhaHi - mAH2 - niA.Hs - itiaHq, 



Gl 



{niAi - niA) H5 - roA(l - v ■ v')He. 



We can now use these redefined tensor form factors to write Eqs. 2 and 3 in a simplified form. 
For transitions to states with J = 1/2, these two equations become 



i?f = u{pa,sa) 



7'' Ai + + L42 + B275 + W'^ ^3 + B375 



u{pa„saJ, 



H^ = u(pa,sa) 



^^{Di + EijA +V^'{D2 + E2J5] +v'^'{D3+E3J5 



w(pa,,saJ- 



For transitions to states with J = 3/2, 



Hi = Ua{pA,SA) 



7^L4i+Bi75 A2 + 5275 + w''^ A3 + ^375 



u{pa„saJ, 



H2 = Ua{pA,SA) 



7^hDi +-E175 +v''(D2 + E2j5 + v'f" ( D3 + E3J5 



"(PA6,saJ, 
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PA = iEA,Ph) 

^+ = {Ee, -Pe) 
FIG. 1: Kinematics of tiie dilepton rest frame. 



while for those with J = 5/2, 

+3^^ (^^4 + S475 

For states with natural parity spinors, 



(^j^ (^Ai + Bi75^ + (^A2 + B275 
(^7^ (^Di + Si 75^ + 1-" (^D2 + E2-/5 



■375 



and for states with unnatural parity, 



A^ = 




CgFi, 


B, = 






= 


CioFi, Ei — 




A, = 




CgGi, 




9^ 


CgFi, 



Di — —CioGi, Ei — C'loFi. 
B. Decays Rates and Forward-Backward Asymmetries 



The decay rate is given by 

dr 



(46) 



(47) 



(48) 



(49) 



(50) 



For the case of unpolarized baryons, is the squared amplitude averaged over the initial polarization and summed 
over the final polarizations. 

For dileptonic decays, the squared average amplitude is 



^„ G 



\M\ 



2 _ ^F^em |T/-,,TA*|2 ( TTf^i^ i 



2%2 



(51) 
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For unpolarized leptons, the leptonic tensors are 

Lai^u = ^^P'^^''^^ = 4 [P+I^P-^ + P+^P-,. - [P- ■ P+ + mj)g^,u] , (52) 

spin 

Lh^u = ^^P^^i"^^ = 4 [p+i^P-i' + P+vP-^, - {p- ■ P+ - m})gi,^] , (53) 

spin 
spin 
spin 

The hadronic tensors H^^ are 

= E^i''^^^ = J2^2^^2, (56) 

pol pol 

= Ei^r^^^2^ h^'' = Yh^^h^- (57) 

pol pol 

The most general Lorentz structure for each of these hadronic tensors is 

H'^-' = -a^g''" + /Si+Q^Q-" + pi.Q^" + fSL+q^Q" + pLg'^q" + ij^e^^^^Qc^q^, (58) 

where Q = pa^ + Pa and q = pa^ — Pa is the 4-momentum transfer. 

We carry out our calculations in the dilepton rest frame (see Fig. 1). In this frame 

^A, = — ^(1-r + s), EA = —^{l-r-s), 



"tA5 / TP '"'Ab fi inAb nr—p: 

Ph = ~^y^~' Ei = —^VS, Pi = ^^y's'ipis 

Ph is the 3-momentum of either baryon in this frame. In addition, we have defined 

s = q^/m\^, r = m,\/m,\^, me = me/niA,, 
(t){s) = (l-r)2 -2(1+ r)s + s^ V(s) = 1 - 4m|/s- 
Performing the contractions, the differential decay rate becomes 

cPT TTiA.Glal^ 



dM~z ~ 2i37r5 i^*"^**^'' VmmMSJ), (59) 
where z = cos 6. In these decays, 4m^ <§<{! — yjrf' and — 1 < -2 < 1. The normalized rate Tq{s, z) has the form 

Tq{s, z) = Io{s) + zli{s) + 5^X2(5), (60) 



where 



and 



Io(s) = a" A, + + a>>B^ + + 0i_B+_ + 0'L+B_+ + (61) 



Aa = 4ml^{2mj+s), A++ = 2m\j{s) , 

Be, = 4mi^(s - 6me), B++ = 2m\^ (0(s) + 4m| (2(1 + r) - s)) , 
S+_ = S_+ = %m\jh\{\ - r), B__ = SmX^mjs. (62) 

The terms proportional to z and z^ are 

Ms) = Amij^^{mS)ir + 7% (63) 
X2{s) = -2mi^0(s)^(s)(/3;+ + /3^+), (64) 
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respectively. Integrating out the z dependence from Eq. 59, the decay rate becomes 



ds 



where 



7eo(s)=Io(s)+ 3X2(5). 



(65) 



(66) 



The functions Xq and I2 are given in Eqs. 61 and 64, respectively. The explicit forms of the coefficients a, I3±±, and 
7 are given in C. 

The forward-backward asymmetry (FBA) is defined as 



Afb{s) = 



1 



dV/ds 



Jo dsdz J_i 



dz 



d^T 



dsdz 



Using Eqs. 59 and 60, we obtain 



Afb{s) 



Ms) 



2 [Ms) + il2(s)] ' 



(67) 



(68) 



where I\ is given in Eq. 63. Since this observable is written as a ratio, it might be expected to be less dependent on 
the form factors than the differential decay rate. HQET considerations may make it even less modcl-dcpcndcnt [57], 
so that it may provide a (somewhat) model-independent way of extracting the Wilson coefficients. Since Afb also 
depends on the chirality of the hadronic and leptonic currents, this observable is also sensitive to any new physics 
effects beyond the Standard Model. 



III. HQET 

Heavy quark effective theory (HQET) is an effective tool for examining the phenomenology of hadrons containing 
a single heavy quark. It has been applied to hadronic matrix elements in many processes at higher and higher order 
in the l/mg expansion, with mg being the mass of the heavy quark. In this section we present the leading-order 
relationships among the HQET form factors and those presented in Section II A for the heavy to light transitions 
Ab — ^ A^*) . The properties of the spinor-tensors that are used to represent the daughter baryons have been discussed 
near the start of section II A. 

To leading order in HQET, transitions between a heavy baryon and a light one are described in terms of only two 
form factors. For any current operator F, the hadronic matrix elements involving tensor states can be written as [86] 



(A(*)(pa(.,) I sTb I K,{v))=u^,...^^{pf,(,,)M^^-^-Tu{v) 



(69) 



where v is the velocity of the A^ baryon and M'^^- -^" is the most general tensor consistent with HQET that can be 
constructed from the available kinematic variables. We may not use any factors of 7^' , , or t^^^'^j in constructing 
^Mi - Mn. therefore, it must take the form 



(70) 



For the case of a pseudo-tensor daughter baryon, the matrix element has the same form as Eq. 69 except that M^^ - ''" 
is now a pseudo-tensor. A pseudo-tensor when sandwiched between spinors can be constructed by multiplying an 
ordinary tensor by 75; thus, for transitions involving a pseudo-tensor spinor. 



c'r\v-p^,,,)+^ct\vPM')) 



75- 



(71) 



Using these forms of the matrix elements, the form factors defined in Section II A are found to satisfy the following 
relations in the limit when the 6 quark is infinitely heavy. 
For any state with = 1/2+ , 



H3^H4^ 0, F2^G2 = -H2 = 2i 



(0) 



(0) 



^(0) 
S2 ) 



(72) 
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while for = 1/2 , we have 



For JP = 3/2- 



For JP = 3/2+, 



For JP = 5/2+, 



F3 = G3=H3 = Hi = 0, F2=G2=H2 = -2Cf , 



i^i = - 



, Gi = —Hi = — 



a(0) ^(0) 
Si ~ S2 



-F3 = G3 = H3 = Fi = Gi = Hi = Hq = He =0, ^2 = ^2 = H2 = —2^2^\ 



Gi = -H, 



Si ~ S2 



^3 = G3 — H3 — F4 — G4 — H4 — H5 = He = 0, F2 = G2 = —H2 = 2^2^\ 



(73) 



(74) 



(75) 



(76) 



Again, we stress that the above relationships are for heavy to light transitions only, and are valid in the limit of an 
infinitely heavy b quark. 



IV. THE MODEL 
A. Quark Model 

In the models considered here, a baryon state takes the form 

Sl, 82,85 

X I [Qi{Pi,si)q2{P2,S2)q{p,Sg)]j^J, (77) 

where Pp = -^{pi—P2),Px = "^(Pi +P2 — 2p) are the Jacobi momenta, C"^ is the antisymmetric color wave function, 
J is the total angular momentum, Eb^ is the energy of the baryon Bq, and the notation [?i929]b denotes the flavor 
wave function of the baryon Bq. For example, the flavor wave function for Ag has the form 

[udq]f^ ^-^{ud—du)q, (78) 

which is antisymmetric in the quarks u and d. When coupled with the appropriate Clebsch-Gordan coefficients, the 
ket I [• • -js^) becomes the flavor-spin wave function for the baryon Bq. x is the momentum space wave function, and 
is given by 

XM1M2M3 ' ^1 Ml i M2 ^s'p Mo - Mi MLSiMs''^I^iMLinplpnxl\)i {Pp,Px), (79) 

i 

where is the total spin of the pair of spectator quarks, Mn is its projection. Si is the total spin of the quarks in 
the baryon, and Ms is its projection. In addition, 

i>LMLn,lpnxlxiPp^P\) = <^it™pi^™x "^Vp^p ("p! ^p)'?^"A;^m^ ("A! Pa)- (80) 

The r]f'' are expansion coefficients that are determined by a diagonalization of a quark model Hamiltonian and 
the C''^"mij2m2 ^-"^^ Clebsch-Gordan coefficients. The combination of momentum, spin and flavor wave functions is 
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symmetric. The functions 4>nlm 

{a;p) are the basis functions. The basis used in this work is the harmonic oscillator 

basis, 



<l>nim{a;p) = exp (-pV2q') 4+'/' (pV«') yim{p), 
where a is the length parameter, which is determined by a minimization of the quark model Hamiltonian, 



(»)'(-!)" 

Q;i+3/2 



2r{n + Ij 
T{n + 1 + 3/2) 



1/2 



(81) 



(82) 



are the associated Lagucrrc polynomials and are the solid harmonics. 
The states described by Eqs. 77-81 are obtained from a variational diagonalization of a quark model Hamiltonian. 
The Hamiltonian used in this model takes the form [87] 



Ki is the nonrelativistic kinetic energy of the ith quark, 



Ki = rrii 



A. 

2171,; 



The spin-independent part of the potential has the linear-plus-Coulomb terms, 



= Ihr - - 



(83) 



(84) 



(85) 



where rij = \fi — rj\. The spin-dependent parts describe the hyperfine and spin-orbit interactions. The hyperfine 
piece, which contains contact and tensor terms, is written as 



ate 



3{Si ■ fij){Sj ■ fij) 



Si • Sj 



and the ad-hoc spin-orbit potential is chosen to be 



m| (p2 + A2) 



L-S. 



(86) 



(87) 



Here ruB^ is the sum of the masses of the three quarks that make up the baryon Bq, p and A are Jacobi coordinates 
conjugate to Pp and p\, respectively, L is the total orbital angular momentum, and S is the total spin of the baryon. 
The parameters Cqqq, b, Q-coul-, ctcon, cttens, ctso, and rrii are obtained from a fit to the experimental spectrum of 
baryon states. 



B. Form Factors 

In order to extract the form factors, we need to calculate matrix elements that take the form 
{A{pA,s')\sTb\Ab{0,s)) = 3-3/2 ^/2^y2;^ j d^p'pd''p'^Sppd^pxC^*C^ 

X E E xiiCe^,{i^p,i^x) 



where 



x{q[q'2S I sTb \ qmb)Xs^s26,iPp^Px) 
{q[q'2S I sTb \ qiq2b) = {q'-^q'^ \ qiq2){s \ sTb \ b). 



(88) 



Rare dileptonic decays of A(, in a quark model 



13 



The matrix element {q'iq2 \ 91 92) gives (5-functions in spin, momentum and flavor in the spectator approximation. 

We calculate the matrix elements of Eq. 88 using two approximations. In the first approximation, we use the 
method outlined by Pervin et al. in [84, 85] . The form factors are extracted by performing a nonrclativistic reduction 
of the quark spinors, keeping terms up to 0(l/mg). The hadronic matrix elements are then calculated analytically 
by using single component wave functions. The matrix elements calculated using the quark model wave functions are 
then used to extract the form factors defined in Section II A. Due to the choice of basis, the resulting form factors 
have the general form of a polynomial times a Gaussian. Once the form factors have been extracted, the polynomials 
are truncated to the zeroth order in the daughter baryon momentum pA. In fact, although the analytic expressions 
for the form factors shown in [84, 85] were shown for single-component wave functions, the results for the decay rates 
shown in those articles were calculated using all components of the wave functions. 

In the second method we use the full multi-component wave functions found from the diagonalization of Eq. 83. 
We also keep the full relativistic form of the quark spinors. Although giving a full analytic treatment of the hadronic 
matrix elements becomes challenging in this case, much of the calculation can still be done analytically and only a 
couple of the integrations need to be performed numerically. The advantage of this approach is that it avoids the 
truncation of the quark currents which may be justified for the h quark but which is much less justifiable for the s 
quark. We will compare the results from the two models with the expectations of HQET. 

We work in the A;, rest frame where the initial quark momenta can be written in terms of the Jacobi momenta as 

Implementing the spectator approximation and integrating over the final Jacobi momenta together lead to 

{k\lTh\K,) = E«^'4^^.i.i'54.,(-1)'^'+'^ 

xf^„;;™, ' "p')W^r;„xUm..; «A'), (90) 



where the coefficients a/j(/j/) are the products of the normalization of the baryon state (the ,J2Eb^ of Eq. 77), the 

expansion coefiicients (the r^" of Eq. 79), and the various Clebsch-Gordan coefiicients that appear in the parent 
(daughter) baryon wave function, and the indices /i(/i') contain all the relevant quantum numbers being summed over 

for the parent (daughter) baryon state. is the spectator overlap. 



This integral can be done analytically and is given in A. W^.^^J'^^^^^' is the interaction overlap, 

x{s{pA +P, v) I sTb I bip,Sq)}(t)n^i^mdl3;P), (92) 

whore /3(') = V^al,'^ is the reduced length parameter for the parent (daughter) baryon, mj\ = m,, + 2rnq. nig is the 
mass of the strange quark and niq is the mass of each light quark. Using the basis functions in Eq. 81, we obtain 

W'r;™,z,m,s; {ax,ax') = J d pexp (^-^ - Cny j yiym.XP ) 

Ix + h ( P" 



x{s{PA+P,Sq') I sTb I b{p,Sg))jCnx ' (^^j yixm^{p), 

where j/ = {2mq/mA)pA + P- The angular dependence in the exponential is eliminated by making the substitutions 

p = k + CPA, = k + dpA, (93) 



where 
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and aw = \/{a'j^ + a'y)/2. This leads to 

nv/vmvV^ ^ _ ( Pa \ /w3. -«ifcV'^' + 5* /^^^ V* (^) 



X {S{PA + P, s,') I sTb I S,))CnV^ (^^ j J^,,™, (p), (95) 

where kx = a\\> /ay ax. The major difference between the two approaches to extracting the form factors is the way 
in which Eq. 95 is handled. 



1. SCA model 

As an example of how the matrix elements arc compiitcd in the SCA approximation, wc examine the ground 
state-to-ground state transition. Using single component wave functions, the spectator overlap is 

C/orK,«P')=(^) ' , (96) 



where Uppi = yia'^ + a^,)/2. The interaction overlap, Eq. 95, becomes 

X j (fke-^"^''\s{pA+p,Sg,) I sTb I b{p,Sg)). (97) 

In this approximation, a nonrelativistic expansion of the spinors in the quark current {s{pA +p,Sq') \ sTb \ b{p,Sg)) 
is carried out. For the resulting matrix element, only terms up to ) are kept. For example, the current 

{s{pa +P, — 1/2) I S7+756 I b{p, +1/2)) is reduced to the form 



{s{pA +P, -1/2) I S7+756 I b{p, +1/2)) = — 

mbrus 

where 



-(pA)+3^n(p) + ^^3^22(p) 



(98) 



7+ = --^(71 +172) , {pA)+ = --^{{pA)i+i{pA)2)- (99) 

The yim{p) can be written in terms of yim{k) by use of an addition theorem. A quark current will thus have the 
general form 

2 e 

{s{pA+P,s,,)\srb\b{p,s,))=J2 E {41iX+br-Zy)yirne{k), (100) 

e=o me=-e 

where ar;eme ^^'^ br;eme depend on mt, rus, and pa; these coefficients are truncated at 0{ _^^_^ ). From Eq. 93, 

p^ = (?Pa + 2cpA ■ k + k"^. The dot product is proportional to X^(— l)'"3^i-m(cpA)3^iTO(fc)- Thus, the integrals in Eq. 
97 have the basic form of an integral over a product of spherical harmonics times an integral of a Gaussian times a 
polynomial, both of which can be performed analytically. 

The form factors for the vector and axial vector currents have the same forms as those published in [84]. For the 
tensor currents, the form factors we obtain are given in B. 
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2. MCN model 



In the MCN model, we keep both the full quark model wave function and the full relativistic form of the quark 
spinors. We therefore need a more general treatment of the interaction overlap, which is written as (see A) 



(101) 



where the coefficients S^"^^ are given in Eq. A8 and 



x{s{pA+p,s,,) I sTb I b{p,Sg))CnV^ [ ^ ) yx^.,{k). 



(102) 



The full relativistic form of Eq. 98 is 



isip, +p, -1/2) I 1 m +1/2)) = -^T^^^WWWT^TWT^' ^'''^ 



where Eb{Es) is the relativistic energy of the b{s) quark. Thus, the current (s(pa +P, Sg') | sTh \ b{p, Sq)) can still be 
written as a linear combination of solid harmonics with maximum i = 2. However, the coefficients of this expansion 
no longer take as simple a form as in Eq. 100. The most general form of the quark current is 

2 e 

{s{pA+P,Sq,)\sTb\b{p,Sq)) = Y^ E ^^T^/r:Sk,pA)ye^M (104) 

e=o me=-e 

The coefficients contain terms of the type k'^ E^^^^{Eb + mb)'^^^^Et^^^{Es + m^)^-^/^, where iV is a nonnegative 
integer. Both Et, and Eg are functions of pa ■ k, which complicates a full analytic treatment of the matrix elements. 
The normalized Lagucrrc polynomials in Eq. 102 arc also functions of p\ ■ k. The coefficients ^ and the normalized 
Laguerre polynomials can be expanded in terms of spherical harmonics so that the momentum and angular integrals 
can be performed separately, with the angular term being done analytically and the momentum integrals being treated 
numerically. The details of this semianalytic treatment of the matrix elements are given in A. 



V. RESULTS 



A. Form Factors 



The parameters for the quark model wave functions used in this work are taken from [87] and are given in Tables 
I and II. For ease of use, all of the form factors calculated using the MCN model are parametrized to have the form 

F{s) = {ao + a2pl + a4pX)e^pl-^^], pj, = 1^,/^, (105) 

where pa is the daughter baryon momentum in the Ab rest frame. The parameters ao, a2 and 04 for the vector and 

axial vector form factors arc given in Tabic III. The corresponding parameters for the tensor form factors arc given 
in Table IV. We have also tried a parametrization that included a term in p\, but the coefficient of that term was 
significantly smaller than the coefficient of the p\ term for all form factors. 

Figure 2 shows a comparison of the vector form factors obtained using the two models for transitions to the ground 
state (Fig. 2(a)), the ffist radial excitation (Fig. 2(b)), the A(1405) (Fig. 2(c)) and the A(1820) (Fig. 2(d)). In the 
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TABLE I: Hamiltonian parameters obtained from a fit to a selection of known baryons. 



(GeV) 


rris 
(GcV) 


rric 
(GcV) 


nib 
(GcV) 


b 

(GcV") 






aso 
(GcV) 




Cqqq 

(GcV) 


0.2848 


0.5553 


1.8182 


5.2019 


0.1540 


^ 0.0 


1.0844 


0.9321 


-0.2230 


-1.4204 



TABLE II: Baryon masses and wave function size parameters, Op and a\, for states with different J^. All values are in GeV. 



State, 


Experiment 


Model 


ax 


Qp 


A6(5620) 1/2+ 


5.62 


5.61 


0.443 


0.385 


A(1115) 1/2+ 


1.12 


1.10 


0.387 


0.372 


A(1600) 1/2+ 


1.60 


1.71 


0.387 


0.372 


A(1405) 1/2" 


1.41 


1.48 


0.333 


0.320 


A(1520) 3/2" 


1.52 


1.53 


0.333 


0.308 


A(1890) 3/2+ 


1.89 


1.81 


0.325 


0.303 


A(1820) 5/2+ 


1.82 


1.81 


0.325 


0.303 



SCA model, all of the kinematic dependence of the form factors takes the form of a Gaussian in the momentum of 
the daughter baryon, calculated in the rest frame of the parent. Thus, all of the SCA form factors have magnitudes 
that monotonically increase as s is increased. In contrast with this, form factors obtained in the MCN have additional 
dependence on the momentum of the daughter baryon, as parametrized in Eq. 105. Depending on the relative sizes 
of oo and 02 in this expression, the shapes of the MCN form factors can be quite different from those obtained in the 
SCA model. The most striking differences occur for the transitions to the radially excited state, shown in Fig. 2(b). 

For all transitions examined, Fa (and F^, where appropriate) is very small in both models, consistent with the 
expectations of HQET. For states with natural parity, Fi is the largest form factor, while for states with unnatural 
parity, F2 is largest. Among the axial form factors, Gi and G2 are the largest, independent of the transition considered, 
and G3 and G4 (where appropriate) are again much smaller. This can be seen by examining the coefficients given 
in Table 111. The SCA form factors all follow a similar trend. Among the tensor form factors. Hi and H2 are the 
largest, with the other form factors all being significantly smaller, in both models, for all transitions considered. 



TABLE III: CoeflScients in the parametrization of the vector and axial-vector form factors obtained in the MCN approach. 





a„(GeV-") 


Fi 


F2 


Fs 


Fi 


Gi 


G2 


G3 


G4 




fflo 


1.21 


-0.202 


-0.0615 




0.927 


-0.236 


0.0756 




Afe ^ A(1115) 


a2 


0.319 


-0.219 


0.00102 




0.104 


-0.233 


0.0195 






a4 


-0.0177 


0.0103 


-0.00139 




-0.00553 


0.0110 


-0.00115 






ao 


0.467 


-0.381 


0.0501 




0.114 


-0.394 


-0.0433 




A(, A(1600) 


a2 


0.615 


-0.281 


-0.0295 




0.300 


-0.307 


0.0478 






04 


0.0568 


-0.0399 


-0.00163 




0.0206 


-0.0445 


0.00566 






ao 


0.246 


-0.984 


0.118 




1.15 


-0.874 


0.00871 




Ab A(1405) 


02 


0.238 


-0.0257 


0.0237 




0.260 


-0.0264 


-0.0196 






a4 


0.00976 


0.0173 


-0.000692 




-0.00303 


0.0159 


-0.000977 






ao 


-1.66 


0.544 


0.126 


-0.0330 


-0.964 


0.625 


-0.183 


0.0530 


Ab A(1520) 


a2 


-0.295 


0.194 


0.00799 


-0.00977 


-0.100 


0.219 


-0.0380 


0.0161 




04 


0.00924 


-0.00420 


-0.000365 


0.00211 


0.00264 


-0.00508 


0.00351 


-0.00221 




ao 


-0.460 


1.33 


-0.232 


0.0485 


-1.71 


1.14 


0.0193 


-0.0153 


Ab A(1890) 


a2 


-0.271 


0.00439 


-0.0315 


0.0140 


-0.284 


0.00990 


0.0374 


-0.00770 




04 


-0.0116 


-0.0149 


0.000345 


-0.00218 


-0.00146 


-0.0134 


-0.000343 


-0.000236 




ao 


2.48 


-0.952 


-0.202 


0.0810 


1.25 


-1.12 


0.355 


-0.143 


Ab ^ A(1820) 


a2 


0.362 


-0.238 


-0.0119 


0.00573 


0.122 


-0.272 


0.0446 


-0.0197 




a4 


-0.00639 


0.00224 


0.000303 


-0.000169 


-0.00134 


0.00303 


-0.00103 


0.000440 
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TABLE IV: Coefficients in tiie parametrization of the tensor form factors obtained in the MCN approach. 





a„(GeV ) 


Hi 


H2 


Hz 


Hi 


Ht, 


Ha 




fflo 


0.936 


0.227 


-0.0757 


-0.0174 






At A(1115) 


02 


0.0722 


0.265 


-0.0195 


-0.00986 








at 


—0.00643 


—0.0101 


0.00116 


—0.000524 








ao 


0.121 


0.389 


0.0421 


0.00676 






Ab A(1600) 


02 


0.313 


0.295 


-0.0479 


-0.0242 








04 


0.0101 


0.0550 


-0.00565 


—0.00404 








fflo 


-1.13 


-0.872 


0.00645 


-0.112 






At A(1405) 




-0.256 


—0.0241 


-0.0197 


-0.00215 










0.00288 


0.0158 


—0.000965 


0.00151 








fflo 


-1.08 


-0.507 


0.187 


0.0772 


-0.0517 


0.0206 


At A(1520) 


a2 


-0.0732 


-0.246 


0.0295 


0.0267 


-0.0173 


0.00679 




04 


0.00464 


0.00309 


-0.00107 


-0.00217 


0.00259 


-0.000220 




fflo 


1.68 


1.13 


0.0214 


0.198 


—0.0147 


0.0331 


At A(1890) 




0.280 


0.00710 


0.0380 


-0.00103 


-0.00818 


0.00674 




04 


0.00154 


-0.0134 


-0.000450 


-0.00155 


-0.000234 


-0.00239 




fflo 


1.55 


0.830 


-0.355 


-0.160 


0.143 


-0.0581 


At ^ A(1820) 




0.0959 


0.298 


-0.0446 


-0.0327 


0.0198 


-0.0205 




04 


-0.00427 


-0.0000926 


0.00103 


0.000739 


-0.000441 


0.00221 




(c) 



(d) 



0.4 
0.2 

=0 
o 

0-0.2 
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t 
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FIG. 2; Comparison of vector form factors calculated using the SCA and MCN models. These form factors are for decays to 
(a) A(1115), = 1/2+, (b) A(1600), = 1/2+, (c) A(1405), = 1/2", and (d) A(1890), = 3/2+. 
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1. Comparison with HQET 

It is of considerable interest to compare the two sets of form factors we obtain with the expectations of HQET. In 
the following paragraphs, we carry out this comparison. 

1.1. jP^l/2+ 

For decays to daughter baryons with = 1/2"'', the leading order HQET predictions are shown in Eq. 72. Figures 
3(a) and (b) show the form factors for decays to the ground state obtained using both the SCA and MCN models, 
respectively. From the figures, we see that both models satisfy the leading order predictions of HQET. In both models, 
F2 ~ G2 through most of the kinematic range. We also note that Gi and Hi are virtually indistinguishable in both 
models. 

The expressions for the leading order HQET predictions can be inverted to give 

= Fi+F2/2 = Gi-G2/2^Hi+H2/2, 
C^") = F2/2 = G2/2 = -H2/2. 

It is useful to extract the independently from the vector, axial vector and tensor form factors, and compare the 
three different forms obtained in this way. Figure 4 shows the three extractions for these form factors from both the 
SCA and MCN models for transitions to the ground state. As can be seen, the form factors obtained using the axial 
vector and tensor form factors are virtually identical in both models. The set obtained from the vector form factors is 
also in very good agreement with the other sets as well. This shows clearly that the form factors obtained from both 
the SCA and MCN models satisfy the leading order expectations of HQET. However, the curves for the ^'"^ from 
the three extractions should not be expected to be identical, since the expressions for the Fi, Gi and Hi in terms of 
universal HQET 'Isgur-Wise' type functions will receive corrections due to the finite mass of the b quark. This holds 
for all of the angular momentum states we consider. 

One HQET expectation is that the form factors for A;, — >■ A should be the same as those for Ac — > A, up to terms 
of order l/m^, where rUg is the mass of the heavy quark. At the 'nonrecoil point', or at maximum q^, the CLEO 

collaboration has extracted the ratio £,2^^ /^i''' ~ —0.25 ± 0.14 ± 0.08 [88]. From our analysis, we find that 

^y(o)^^y(o) ^ _o.l66(-0.092), ^^'°V^f = -0.193(-0.113), 
^T(o)^^T(o) ^ _o.l88(-0.108), 

within the SCA (MCN) model, at the nonrecoil point. For both the SCA and MCN models, our values for this ratio 



is smaller in the MCN model than in 



are consistent with the value reported by CLEO. We note here that 

the SCA model for all three (vector, axial vector, tensor) form factor scenarios. 

The HQET form factors for transitions to the first radial excitation obtained from the two models we use are shown 

in Figure 5. For each model, the three possible extractions of the agree with each other reasonably well over the 

entire kinematically allowed range. Not surprisingly, the extracted from the two models are very different. 
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FIG. 4: HQET Form Factors for At ->■ A(1115), = 1/2+ as a function of s = q'^/ml^. The graphs show the HQET form 
factors calculated using the vector, axial vector, and tensor form factors using (a) SCA and (b) MCN models. 




FIG. 5: Same as Fig. 4, but for Af, A(1600), = 1/2 



The shape of the function ^{ in the MCN model can be understood by taking the limit in which the strange 
quark is also heavy. In that limit, the transitions to states with = 1/2+ are described by a single form factor. 
Furthermore, for transitions to the ground state, that form factor, the Isgur-Wise function, is normalized to unity at 
the nonrecoil point. Since the wave function of the radially excited state must be orthogonal to that of the ground 
state, the corresponding Isgur-Wise function for transitions to the radially excited state must vanish at the nonrecoil 
point, up to corrections proportional to the inverse of the heavy quark masses. In the SCA model, the truncated 
form factors still exhibit the expected behavior of being 0{l/mq) near the non-recoil point, but the more precise form 
factors of the MCN model show this behavior more clearly. 

The ratios Ca"'-"/?!" ' at zero recoil for this state for the three extractions within the SCA (MCN) model are 

^no')/^v(o') ^ _o.615(-0.689), ^^^"'V^f = -0.611(-0.633), 
^T(o')/^T(o') ^ _o.590(-0.617). 



For the three form factor scenarios, V'^i" "* 



1.2. 



is smaller in the SCA model than in the MCN model. 



For decays to states with = 1/2 , the leading-order HQET predictions are shown in Eq. 73. These relations can 
be inverted to give 

.(0) _ 



Cr = -{Fi-F2l2) = -{Gi+G2/2)=H,-H2/2, 

^2 



Cf' = -F2/2 = -G2/2 = -i/2/2. 



Figure 6 shows the three versions of the extracted C,^^^ for transitions to this state, for each model. There is good 
agreement among the three sets of form factors in both models. We note here that although Ci'^"'' and Cf"^^ are 
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FIG. 6: Same as Fig. 4, but for Ab -> A(1405), ^ = 1/2" 
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FIG. 7: Same as Fig. 4, but for Ab A(1520), = 3/2" 



indistinguishable over the entire kinematic range in the SCA model, they begin to deviate slightly from each other as 

C^^°^ and 



they approach zero recoil in the MCN model. In both models, (2^^"^ ^^'^ Cz^^^ virtually identical. The values for 



the ratio C2/C1 zero recoil for transitions to this state are 

^y(o)/^v^(o) ^ _o.555(-0.667), 
^T(o)^^T(o) ^ _o.517(-0.626). 



within the SCA (MCN) model. From the above, we see that C2°VCi°' 
scenarios. 



-0.499(-0.610), 



is larger in the MCN model for all three 



1.3. J' 



3/2- 



The leading-order HQET predictions for decays to states with = 3/2 arc shown Eq. 74. By inverting these 
expressions, we find that 

^ Fi+ F2/2 = Gi - G2/2 = Hi+ H2/2, 
^i'^ = F2/2 - G2/2 = -H2/2. 

The three sets of HQET form factors for transitions to this state for each model are shown in Figure 7. The three 
extractions agree with each other throughout most of the kinematic range, in both models. For ^i^^ /^[^^ at the 
nonrccoil point, we find 

^ni)/^vii) ^ _o.l39(-0.196), e^'^V^f = -0.172(-0.245), 
^T(i)^^T(i) ^ _o.lll(-0.19G). 

We see that the magnitude of this ratio is larger in the MCN model for all three form factor scenarios. 
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FIG. 8: Same as Fig. 4, but for At A(1890), = 3/2+, 

(a) (b) 





FIG. 9: Same as Fig. 4, but for Af, A(1820), r = 5/2^ 



1.4. J' 



3/2~ 



For decays to states with = 3/2+, the leading-order HQET predictions are given in Eq. 75. Inverting these 
expressions leads to 

Cf^ = -{Fi-F2/2) = -{G, + G2/2)^H,~H2/2, 



/■(I) 

^2 



-F2/2 = -G2/2 = -H2/2. 



Figure 8 shows the three extractions for the (^^'^'^ for transitions to this state obtained in both models. There is good 
agreement among the three sets of form factors for both models. In the SCA model, the three sets of form factors are 
nearly the same over much of the kinematic range. However, in the MCN model we find that the cj^^ for all three 
scenarios are almost indistinguishable over the entire kinematic range, but C2''^^ deviates from the other two as they 
approach the nonrecoil point. The ratio C2^VCi^'' ^-t zero recoil for the three form factor scenarios are 

^v{i)/^v{i) ^ _o.353(-0.591), Ca'^^'Vcf = -0.285(-0.501), 
^T(i)^^T(i) ^ -0.293(-0.509), 

within the SCA (MCN) model. The magnitude of this ratio is larger in the MCN model for all three scenarios. 
1.5. jP = 5/2+ 

For decays to daughter baryons with = 5/2+, the leading-order HQET predictions are shown in Eq. 76. These 
expressions can be inverted to give 



(2) 



-I- F2/2 = Gi - G2/2 = Hi+ H2/2, 



^^'^ = F2/2^G2/2^-H2/2. 
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TABLE V: SM and SUSY values for the Wilson coefficients. In tlie SUSY model we use, only C7, Cg, and Cio get modified 
from their SM values. 





Ci 


C2 


C3 


Ci 


a, 


Ce 


C7 


Ca 


Clo 


SM 


-0.243 


1.105 


0.011 


-0.025 


0.007 


-0.031 


-0.312 


4.193 


-4.578 


SUSY 














0.376 


4.767 


-3.735 



The three sets of HQET form factors for transitions to this state are shown in Figure 9. In both models, the three 
sets of form factors are nearly identical at small s. In the SCA model, ^ ' and ^ ' remain indistinguishable over 
the entire kinematic range while deviates as the nonrecoil point is approached. It can also be seen that ^^^^^ 

and ^2 remain fairly close over the kinematic range with ^2 deviating from the others as they approach zero 
recoil. In the MCN model, the sets of form factors obtained from the vector and tensor form factors remain nearly 
the same throughout the entire kinematic range with those obtained from the axial vector form factors deviating as 
they approach zero recoil. The ratio ^2^VCi^' nonrecoil point for the three scenarios are 



-0.130(-0.238), C^^^Vcf^^^ 



-0.172(-0.309), 



^T(2) 



7^p) = -0.056(-0.211), 
for the SCA (MCN) model. The magnitude of this ratio is larger in the MCN model for all three scenarios. 



B. Dileptonic Decays 

In this section, we present differential decays rates, branching ratios (BRs), and FBAs for dileptonic decays using 
the two sets of form factors that we have extracted. We present results obtained using Wilson coefficients that have 
been calculated in the standard model (SM) , with both sets of form factors. We also examine one scenario that arises 
beyond the SM, namely a supcrsymmctric (SUSY) extension to the SM, but there we use the MCN form factors 
exclusively. In our numerical calculations, the SM values of the Wilson coefficients are taken from [55] and the SUSY 
values are taken from [53]. These values are presented in Table V. In [55], the Wilson coefficients are evaluated using a 
naive dimensional regularization scheme at the scale = 5.0 GeV. The top quark mass is taken to be rrit = 174 GeV 
and the cut-off A;^^ = 225 MeV, where MS denotes the modified minimal subtraction scheme. The SUSY model used 
here is referred to as SUSYI in Ref. [53] . SUSYI corresponds to the regions of parameter space where supersymmetry 
can destructively contribute and can change the sign of C7, but contributions from neutral Higgs bosons are neglected. 

Since the rates and FBAs for the e and ji channels are essentially the same, in what follows, we present the results 
for the /i and r channels only. 



1. Decay Rates 

The branching ratios predicted for the /z and r channels are presented in Tables VI and VII, respectively. Each 
table displays the results for the SM calculations using two models for the form factors, as well as one SUSY scenario 
with the MCN form factors. In addition, results obtained omitting and including the long distance (LD) contributions 
are presented. For ease of discussion, we will refer to the results obtained in the SM as SMI for the SCA form factors, 
and SM2 for the MCN form factors. SMla and SM2a will refer to results with LD contributions omitted, while SMlb 
and SM2b will refer to results with LD contributions included. Finally, SUSYa and SUSYb will refer to the results 
obtained using Wilson coefficients from the supcrsymmctric extension to the standard model discussed above, with 
the SUSYa (SUSYb) results obtained when LD contributions are omitted (included). We also compare our model 
predictions with those made by other authors using SM Wilson coefficients, within the framework of light cone sum 
rules (LCSR) [53, 54], QCD sum rules (QCDSR) [56], and a multipole model (PM) [56], as well as with the recent 
experimental results from the CDF Collaboration [82]. Graphs of the differential decay rates dT/ds are shown in Figs. 
10-15. 

For decays to the ground state in the fj, channel, the branching ratios we obtain in the SM scenario are smaller 
than those obtained by other authors [53, 54, 56]. In the t channel, our results in models SMla and SM2a are smaller 
than the LCSR predictions for the ground state, but are comparable to the results obtained using QCDSR and PM. 
However, the predictions of the SMlb and SM2b models are much smaller than those of other models in this channel. 
The predictions of SM2 are larger than SMI for decays to all A^*^ for both channels. Note, too, that SMI predicts 
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TABLE VI: Branching ratios for At H> A^*' /i"*"/^ in units of 10 ^. The numbers in the column labeled SMI are obtained using 
the SCA form factors with standard model Wilson Coefficients. The numbers in the column labeled SM2 are also obtained 
using SM Wilson coefficients, but the MCN form factors. The numbers in the column labeled SUSY are obtained using the 
MCN form factors with Wilson coefficients from a supersymmetric scenario. The column labeled LD refers to the long distance 
contributions of the charmonium resonances, with 'a' indicating that these contributions have been neglected, and 'b' indicating 
that they have been included. In this table, it is assumed that the A(1600) is the first radial excitation. The lifetime of the A^ 
is taken from the Particle Data Listings [89]. 



btate, J 


LD 


SMI 


aM2 


SUSY 


Aslam et al. [53J 


Wang et al. [54] 


Chen et al. [56] 
QCDSR PM 


TT^ « I i Tool 

tixperiment [82] 


A(1115) 1/2^ 


a 


0.60 


0.70 


1.0 


5.9 


6.1 


2.1 1.2 


1.73 ± 0.42 ± 0.55. 




b 


21 


32 


32 


39 


46 


53 36 




A(1600) 1/2+ 


a 


0.027 


0.32 


0.53 












b 


2.6 


35 


35 










A(1405) 1/2" 


a 


0.094 


0.21 


0.32 












b 


5.9 


19 


19 










A(1520) 3/2^ 


a 


0.13 


0.21 


0.34 












b 


14 


24 


24 










A(1890) 3/2+ 


a 


0.018 


0.097 


0.17 












b 


1.3 


5.8 


5.9 










A(1820) 5/2+ 


a 


0.013 


0.082 


0.15 












b 


0.84 


4.6 


4.7 











TABLE VII: Branching ratios for A^ — ;> A'^*V+r in units of 10 The columns are labeled as in Table VI. 



State, J-^ 


LD 


SMI 


SM2 


SUSY 


Aslam et al. [53] 


Wang et al. [54] 


Chen et al. [56] 
QCDSR PM 


A(1115) 1/2+ 


a 


0.22 


0.22 


0.38 


2.1 


2.4 


0.18 0.26 




b 


0.59 


0.68 


0.86 


4.0 


4.3 


11 9.0 


A(1600) 1/2+ 


a 


< 0.01 


< 0.01 


0.016 










b 


0.033 


0.12 


0.13 








A(1405) 1/2" 


a 


0.023 


0.030 


0.055 










b 


0.12 


0.22 


0.25 








A(1520) 3/2" 


a 


0.013 


0.016 


0.031 










b 


O.ii 


0.18 


0.20 








A(1890) 3/2+ 


a 


< 0.01 


< 0.01 


< 0.01 










b 


< 0.01 


< 0.01 


< 0.01 








A(1820) 5/2+ 


a 


< 0.01 


< 0.01 


< 0.01 










b 


< 0.01 


< 0.01 


< 0.01 









that decays to the ground state dominate the rare decays of the Af,, but SM2b indicates that decays to the first radial 
excitation are the dominant mode, with the decay rate to the A(1520) being similar in magnitude to the decay rate to 
the ground state. In addition, the decay rate to the A(1405) is only slightly smaller than the decay rate to the ground 
state. These results imply that searches for rare decays of the A;,, such as those planned by the LHCb collaboration, 
should include excited final states, as the decay rates to these states can be sizable. We note that this is consistent 
with the current experimental status in the rare decays of the B meson, where decays to the K and K* account for 
less than half of the inclusive dileptonic decay rate. In both the fj, and r channel, the SUSY BRs are larger than the 
SM results without LD contributions. 

J'' = 1/2+ 

Figs. 10(a) and (b) show the differential decay rates for decays to the ground state (A(1115)) in the fx and r channels, 
respectively. In these graphs, the solid curves are obtained from SMI, the dashed curves are from SM2, while the 
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FIG. 10: dV/ds for (a) A;, — >■ A(1115)/i+/:i" and (b) Kb A(1115)r"'"r~ without and with long distance (LD) contributions 
from charmonium resonances. The sohd curves represent rates obtained from SMI, the dashed curves are from SM2, and the 
dot-dashed curves are from SUSY. 
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FIG. 11: Same as Fig. 10 but for At A(1600). 



dot-dashed curves are from SUSY. The resonance contributions enhance the rates for both decay channels, but the 
enhancement is much larger in the dimuon case. The SMla and SM2a branching ratios are somewhat smaller than 
the experimental measurement, but the SMlb and SM2b branching ratios are significantly larger. Experimentally, 
regions in s around the first two charmonium resonances are vetoed. This is because the non-leptonic decay rates of 
the A;, to A along with a vector charmonium, followed by the dileptonic decays of the charmonium, are much larger 
than the rare decays of interest. Such decays can be described by tree-level diagrams, and thus have much larger 
decay rates than rare decays. To compare our numbers with experiment, wc should follow a similar procedure of 
vetoing the regions around the two lightest charmonium resonances. When this is done, the branching ratios that 
result are essentially identical to those that we report in the 'a' scenarios. Thus, SMla and SM2a are the numbers 
that should be compared with experiment. 

Figure 11 shows the differential decay rates for transitions to the first radial excitation, A(1600). The SMI predic- 
tions shown in Table VI are much smaller than the SM2 predictions for both decay channels. Furthermore, SM2b 
predicts that decays to this state are the dominant rare decay mode of the A;,. The truncations of the quark currents 
and the form factors in SMI have lead to significant underestimates of the rates for decays to this state in both the 
/U and T channels. 



= 1/2- 

Figs. 12(a) and (b) show the differential decay rates to the lowest-lying l/2~ state, the A(1405), assuming that 
it is a three-quark state (there is at least one other suggestion for the structure of this state in the literature [90]). 
From the graphs, it can be seen that SM2 predicts a larger rate than SMI, and this is seen in Tables VI and VII. The 
SUSY rates are significantly larger in both decay channels over the entire kinematic range. 
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FIG. 12: Same as Fig. 10 but for -> A(1405). 
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FIG. 13: Same as Fig. 10 but for A^ A(1520). 

JP = 3/2- 

Figs. 13(a) and (b) show the differential decay rates for decays to A(1520). The curves indicate that SMI and SM2 
make similar predictions for this channel, and this is borne out by the numbers in Tables VI and VII. The SUSY 
rates arc larger in both the /i and r channels. For this state, SM2b predicts that the decay rate into this state is very 
similar to the decay rate into the ground state. 



JP = 3/2+ 

Figs. 14(a) and (b) show the differential decay rates for decays to A(1890). In both the /i and r channels, the 
predicted SM2 rates significantly larger than the SMI rates. The SUSY rates are larger in both decay channels over 
the entire kinematic range. The SM2b rate for the /i channel into this mode is quite a bit smaller than the rate into 
the ground state. Nevertheless, this decay rate is not negligible. 



= 5/2+ 

Figs. 15(a) and (b) show the differential decay rates for decays to A(1820). As with many of the cases discussed, 
SM2 predicts larger rates for decays to this state than SMI, in the /x channel. In the r channel, the predicted rates 
are negligible. The SUSY rates are larger in both decay channels over the entire kinematic range. 



2. Forward-Backward Asymmetries 

The differential forward-backward asymmetries (FBAs) Afb{s) are shown in Figs. 16-21. The key to the curves 
is the same as the differential decay rates. In addition to the differential asymmetries, the zeroes in the asymmetries 
also contain information on the Wilson coefficients, and are therefore of interest. We will discuss these in some detail. 



Rare dileptonic decays of A(, in a quark model 



26 




FIG. 14: Same as Fig. 10 but for A(1890). 
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FIG. 15: Same as Fig. 10 but for Kb A(1820). 



Furthermore, it is useful to introduce the integrated forward-backward asymmetry {Afb) in order to characterize the 
typical value of the FBA. This integrated FBA is defined as 



{Afb) = I AFB{s)ds. 

Jirhl 



(106) 



The integrated FBAs we obtain arc shown in Tables VIII and IX. The column labels have the same meaning as 
with the branching ratios shown in Tables VI and VII. We also compare our results with those of Wang et. al. [54] 



TABLE VIII: Integrated forward-backward asymmetry for Aj — >■ A'*'/u,'''/x . The columns are labeled as in Table VI. 



State, 


LD 


SMI 


SM2 


SUSY 


Wang et al. [54] 


Chen et al. [57] 


A(1115) 1/2+ 


a 


-0.1255 


-0.1272 


-0.1800 


-0.0122 


-0.1338 




b 


-0.1128 


-0.1139 


-0.1582 


-0.0099 




A(1600) 1/2+ 


a 


-0.0663 


-0.0889 


-0.1415 








b 


-0.0575 


-0.0793 


-0.1202 






A(1405) 1/2- 


a 


0.0927 


0.1145 


0.1681 








b 


0.0822 


0.1021 


0.1455 






A(1520) 3/2" 


a 


-0.0611 


-0.0825 


-0.1371 








b 


-0.0566 


-0.0747 


-0.1165 






A(1890) 3/2+ 


a 


0.0410 


0.0586 


0.1116 








b 


0.0431 


0.0580 


0.0985 






A(1820) 5/2+ 


a 


-0.0074 


-0.0271 


-0.0946 








b 


-0.0133 


-0.0311 


-0.0846 
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TABLE IX: Integrated forward-backward asymmetry for A;, — > A'*V"'"r~. The columns are labeled as in Table VI. 



State, J-^ 


LD 


SMI 


SM2 


SUSY 


Wang et al. [54] 


Chen et al. [57] 


A(1115) 1/2+ 


a 


-0.0342 


-0.0330 


-0.0274 


-0.0067 


-0.0399 




b 


-0.0319 


-0.0307 


-0.0259 


-0.0062 


— 


A(1600) 1/2+ 


a 


-0.0119 


-0.0116 


-0.0098 


— 


— 




b 


-0.0106 


-0.0103 


-0.0088 


— 


— 


A(1405) 1/2- 


a 


0.0221 


0.0243 


0.0197 


— 


— 




b 


0.0203 


0.0222 


0.0183 


— 


— 


A(1520) 3/2^ 


a 


-0.0072 


-0.0098 


-0.0086 








b 


-0.0063 


-0.0086 


-0.0077 






A(1890) 3/2+ 


a 


0.0003 


0.0013 


0.0013 








b 


0.0002 


0.0010 


0.0010 






A(1820) 5/2+ 


a 


0.0021 


0.0021 


0.0003 








b 


0.0020 


0.0021 


0.0003 







TABLE X: Zeroes of the forward-backward asymmetry for A;, — > A'*-'/i+/i . The columns are labeled as in Table VI. 



State, 


LD 


SMI 


SM2 


SUSY 


Chen et al. [58] 


A(1115) 1/2+ 


a 


0.106 


0.098 




0.109 




b 


0.088 


0.080 




0.098 


A(1600) 1/2+ 


a 


0.143 


0.102 








b 


0.121 


0.085 






A(1405) 1/2" 


a 


0.132 


0.107 








b 


0.111 


0.088 






A(1520) 3/2" 


a 


0.117 0.522 


0.101 0.525 


0.527 






b 


0.097 0.522 


0.085 0.525 


0.527 




A(1890) 3/2+ 


a 


0.116 0.431 


0.103 0.437 


0.439 






b 


0.102 0.430 


0.094 0.430 


0.439 




A(1820) 5/2+ 


a 


0.138 0.415 


0.112 0.424 


0.438 






b 


0.121 0.422 


0.099 0.430 


0.441 





(LCSR) and Chen et al. [57] (QCDSR). Tables X and XI show the locations of the zeroes in the FBAs. 

J'' = 1/2+ 

Figs. 16(a) and (b) show the differential FBAs for decays to the ground state, in the /it and r channels, respectively. 
The values of the integrated asymmetries, {AFB)^ obtained in the SMla and SM2a models, agree well with the 
values reported by Chen et al. [57] for both channels, but is roughly an order of magnitude larger than the value 
reported by Wang et al. [54]. Our values from the SMlb and SM2b models are also an order of magnitude larger 
than the corresponding value reported by Wang et al. It is very interesting to note that for this final state, the 



TABLE XI: Zeroes of the forward-backward asymmetry for A^ — > A'*V+r . The columns are labeled as in Table VI. 



State, 


LD 


SMI 


SM2 


SUSY 


A(1520) 3/2" 


a 


0.522 


0.525 


0.527 




b 


0.522 


0.525 


0.527 


A(1890) 3/2+ 


a 


0.431 


0.437 


0.439 




b 


0.433 


0.437 


0.438 


A(1820) 5/2+ 


a 


0.415 


0.425 


0.439 




b 


0. 112 


0. 130 


0.110 




integrated asymmetry appears to be largely independent of the model employed. When the results of other authors 
are considered, the results of Chen et al. confirm this conclusion, but the results of Wang et al. do not allow such 
a conclusion to be drawn. The curves of Figs. 16(a) and (b) also indicate that the FBA for decays to this state are 
largely independent of the form factors used. 

The FBAs for decays to A(1600) for the n and r channels are shown in Figs. 17(a) and (b), respectively. In the 
^, channel, the SMI and SM2 models lead to significantly different FBAs over most of the kinematic range, and this 
is reflected in the values of the integrated FBAs. In addition, the locations of the zeroes are quite diff'erent in the 
two models. In the r channel, the curves from the SMI and SM2 models are closer than in the /U channel, and the 
integrated FBAs are essentially identical. 

From Eq. 68, the zeroes of the FBA occur when 

^Fs(so) =O^Ii(so) = 0. (107) 
For decays to states with J = 1/2, the expression for Ii is 

Ti(s) = -l&m\j^/(t){s)i!{s) [Yie{AlEi) + Yic{BlDi)] 
Using Eq. 48, the condition for the position of the zero(s) for decays to states with = 1/2+ is 

So V 2rnAi,FiGi / 

This relation holds for = l/2~ as well. From this relation, we see that so depends on the two combinations of 
Wilson coefficients, 'Re{CjCw) and Re(CgCio), and a ratio of form factors. For this final state, there is a single 
possible zero when LD contributions arc omitted: LD contributions introduce other zeroes, as can be clearly seen in 
Figs. 16(a) and 17(a). Using the form factors obtained in the two models we consider, the values of so are shown in 
Table X for the ji channel. These values are in good agreement with those reported in [58]. 

For decays to states with J = 1/2, apart from the cndpoints and resonance regions, there are no zeroes in the FBAs 
in either channel for the SUSY model we use. The zeroes that occur in the SM scenarios can be traced to the opposite 



(108) 
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signs of C7 and Cg. In the SUSY model that we employ here, C^ and Cg have the same sign, and the condition for 
the zero in Eq. 109 can no longer be satisfied. In addition, in the r channel, even in the SMI and SM2 scenarios, no 
zeros are possible, apart from those induced by the resonance effects. 



= 1/2- 

The FBAs predicted for decays to the A(1405) are shown in Fig. 18. In these figures, we see that the predictions 
from SMI and SM2 are somewhat different, and the locations of the zero are also different (see Table X). As with 
the decay to states with = 1/2+, there are no zeros in the FBA in the SUSY scenario we consider here. There 
are also no zeros, apart from those induced by the resonance contributions, in the r channel for any of the scenarios 
considered. 



S/2- 



The predictions for the FBAs in the decays to the A(1520) are shown in Fig. 19. The SMI and SM2 models give 
slightly different asymmetries in the channel. Even without the LD contributions, the structure of the asymmetry 
arising from these decays is richer than in the decays to states with J = 1/2. 

For states with J = 3/2 



Ii(s) 



is) 



[Re{AlEi) +Re{B*Di)] + 



y-{s) [Re{AlEi) + Re(B|£)i)] + y+{s) [Re{AlEi) + ReiB*Di)] + 



2[Re{AlE4) +Re{B2D4)] 



(110) 




FIG. 20: Same as Fig. 16 but for At A(1890). 



where 0(s) = (1 - r)^ - 2(1 + r)s + and y^{s) = [(1 ± ^/r)'^ - s] /y^. For = 3/2 , the condition for the 
positions of the zeros is 

Re(C9*Cio) = ^Re(C7*Cio) f ' (111) 

where 

X, = -^{F^G,-GjF,)+y-{§o)iF^G4-GjF,) + 

2/+ (so) (f jGi - Gf fi) + 2 (fJG4 - Glf^4) , 

n = -^^F,Gi+y-{so)FiG4 + y+{so)F4Gi+2F4G4, (112) 
r 

As wc can soc in Fig. 19, apart from the resonance region, there are two zeroes for this mode in the /i channel in 
the SM. The zero at the larger value of s is also present in the r channel. This is quite different from the case with 
J = 1/2 where there is only one zero in the fx channel and none for the r. The positions of the zeroes in the /x channel 
are shown in Table X. In the SUSY scenario that we explore, there is only one zero in the fi channel, and it sits at 
s = 0.527 in both SUSY scenarios. In the r channel, there is a single zero in the FBA, and its position is largely 
independent of the model used for the form factors. 



JP = 3/2+ 



Fig. 20(a) shows the /i channel FBAs for the A(1890) mode. There arc some differences between the predictions of 
the SMI and SM2 models, particularly at the larger values of s. These differences appear to be even larger in the r 
channel (Fig. 20(b)), but the asymmetries predicted are quite small in both models. 

For this decay mode, the condition for the zeroes in the asymmetry is the same as in Eq. Ill, but now 

Xr = -^(F[G,-G{F,)+y~{so){FlG,-G^F4) + 
y+(so) {F[G4 - GiFi) + 2 (FJG4 - G^Fi) , 

n = -^FiGi+y-(so)F4Gi+y+(so)FiG4 + 2i=4G4, (113) 
r 

The positions of the zeroes in the FBA for decays to this state are also shown in Table X for the channel and Table 
XI for the r channel. 



= 5/2+ 

In the A(1820) mode, SMI and SM2 again predict slightly different FBAs in both the /x and r channels (see Fig. 
21). In the fj. channel, the SMI prediction for the position of the first zero is noticeably different from the prediction 
from SM2. However, the positions of the second zero are much closer in these two models. 




For J = 5/2, Ii takes the form 



Ms) = -m\jcly^/\s)y/m 

2[Re{AlEi) +Re{BlDi 

Re(A*^4) +Re(5*Li4) 



[Re{A\E{) +Y{j&{BlDi) 



■2[Re(A*^i)+Re(S*D4)] + 



The condition for the zero now becomes 



where, for natural parity, 



Xn = - 



So v2toA(,>2 



{F^Gr - GIF,) - 2 (Ff G4 - GiFi) + 2 {fJG, - GJF,) + 



(114) 



(115) 



FTGa — GJ Fa 



Yi = - 



F\G\ — 2_FiG'4 + 2F4G1 + F4G4. 



(116) 



We conclude this section with two general comments on the zeroes of the FBAs. The zeroes for the FBAs in decays 
to states with spin 1/2, along with the first zero in decays to the states with higher spin that we have considered, all 
lie relatively close to each other, despite the more complicated expression for the location of the zeroes for the states 
with higher spin. However, in Eqs. 112, 113 and 116, terms in F4, G4, Fj and G4 are negligible (they are exactly 
zero in the limit of an infinitely heavy b quark), so that the various expressions for the location of the zeroes all reduce 
to one that is identical to the case of spin 1/2, Eq. 109. Furthermore 



for states with natural parity, or 



FiGi — GJFi 
2mA,FiGi 



F^Gi - GjFi 
2mA,FiGi 



1 + 



"'I 



O 



O 



A, 



QCD 

nib 



Aqcd \ 
rrib J 



(117) 



(118) 



for states with unnatural parity, for the states we have examined. This means that the location of the zero, up to 
corrections O and O ^|^^ or O , is approximately given by 



So 



-2mb 



Re(Cg*Cio) ' 



(119) 



independent of the angular momentum of the final state (at least, up to spin 5/2), and of the form factors. Using 
the SM Wilson coefficients along with the physical mass of the A^ and the accepted mass of the b quark, this gives a 
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0.3 



0.2 



0.1 
J 0.0 
-0.1 

-0.2 



0.00 0.09 0.18 0.27 0.36 0.45 0.54 

s 

FIG. 22: FBA for At — )• A(1520), obtained in model SMla. The solid curve arises from the SCA form factors. The dashed 
curve is obtained by omitting any terms that are 0{l/mb) from the SCA form factors. The dot-dashed curve is obtained by 
allowing F4 to be non-zero. 

value of 0.121. This number, obtained in this simplifying limit, is in surprisingly good agreement with the values of 
the lower zeroes shown in Table X, for all states. Of course, there must be deviations from this simple limit, as the b 
quarks is not infinitely heavy. However, it may be possible to systematically estimate the corrections to the (model 
independent?) value of 0.121. 

The second comment on the zeroes follows immediately from the first: in the limit of an infinitely heavy b quark, 
the location of the zeros is given by Eq. 119 (up to the corrections mentioned), since F4, G4, Fj and G4 all vanish 
explicitly in this limit. This means that in this limit, there can only be one zero in the FBAs if the form factors 
are treated in the strict heavy-quark limit. This also means that the location of this second zero may be sensitively 
dependent on the form factors, the angular momentum of the state being considered and, of course, on the Wilson 
coefficients. Surprisingly, the results shown in Table XI suggest that the most important dependence is the angular 
momentum of the daughter baryon. 

This is illustrated in Fig. 22. In this graph, the solid curve arises from the SCA form factors. The dashed curve 
is obtained by omitting any terms that are ©(l/m^) from the SCA form factors. This means, for instance that _F4, 
G4, FJ and all vanish identically. It's clear that for these truncated form factors, there is only a single zero in 
the FBA. The dot-dashed curve arises when only F^ is allowed to be non-zero, and the second zero reappears in the 
FBA. 



C. Sensitivity to Model Parameters 

One question that can be asked is how sensitive are the results presented above to the choice of parameters in the 
model. The easiest way to assess this is to vary one or more parameters from the values we have used, and note the 
effect on the calculated decay rates, for instance. We have carried out this exercise, changing the mass of the b quark 
by 10% from its value reported in Table I. The result is that the decay rate to the ground state A changes by less 
than 0.02% of the reported value, in either model of the form factors. For decays to other states, the changes in the 
decay rate are similarly small. This suggests that the results that we have obtained are quite robust to changes in 
parameters, within the context of this particular model. It is extremely difficult for us to speculate on how the results 
would change if the model itself were to be changed. We should also emphasize that the parameters reported in Table 
I are obtained in a fit to the baryon spectrum. A change in the mass of the b quark of 10% means that the model 
states are no longer good representations of the observed physical states. 
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VI. CONCLUSIONS AND OUTLOOK 

In this work, we have examined the rare weak dileptonic decays of the A;, baryon using a constituent quark model. 
Analytic and numerical results for the form factors for the decays to groimd and excited states with different quantum 
numbers have been obtained and compared with the leading order predictions of HQET. We have shown that for 
Afe A^*\ the leading order relations among form factors predicted by HQET are satisfied by the form factors for 
both the SCA and MCN models. For transitions considered, MCN model form factors are found to be larger than 
those in the SCA model. 

We have examined the decay rates and FBAs for decays to a number of A final states, in a number of scenarios. We 
have compared our results with the predictions of other authors, as well as the /x-channel measurement reported by 
the CDF collaboration. Our predictions for decays to the ground state are smaller than the experimental value and 
the predictions of other authors in the /i channel. In the r channel, our model results are in agreement with QCD 
sum rules and pole model predictions without LD contributions. However, when LD contributions are considered, our 
results are smaller than those of other authors. We comment on this further near the end of this section. 

The SM2 models predict larger branching ratios for all daughter A decay modes than do the SMI models. SMI 
predicts that decays to the ground state are dominant in the /z channel. In SM2a, this mode is also dominant, but in 
SM2b, decays to A(1600) are dominant. SM2b also predicts that the decay rate to the A(1520) is comparable with 
the decay rate to the ground state, with the decay rate to the A(1405) being only slightly smaller. These results are 
consistent with the current status of the rare dileptonic decays of B mesons, where decays to the two lowest lying 
kaons are insufficient to saturate the inclusive rate, and in fact account for less than 50% of the inclusive rate. These 
results also suggest that it might be prudent for the LHCb collaboration and CDF to search for rare decays in modes 
other than the A(1115). 

There is, however, a very important caveat. Decays in which the of the dileptons is near the mass-squared of the 
two lowest-lying vector charmonium resonances are inaccessible experimentally, as they are embedded in the much 
larger background coming from (tree-level) nonleptonic decays with vector charmonia in the final state. Our results 
obtained omitting the vector charmonia, the SMla and SM2a models, arc therefore closer to experimental reality, and 
these scenarios suggest that decays to the A(1115) are indeed the dominant rare decay mode of the A(,. Nevertheless, 
the SM2a model that uses the more precise calculation of the form factors indicates that there will be a sizable fraction 
of rare decays into excited states of the A. 

For some decay modes, such as decays to the A(1115), the FBAs are largely independent of the model choice for 
much of the kinematic range. For other modes, particularly the A(1600), there are significant differences between the 
predictions of the SMI and SM2 models. The zero that occurs at lower values of s in these FBAs turns out to be 
largely independent of the angular momentum of the daughter baryon, but with some dependence on the form factors. 
Nevertheless, this form factor dependence is surprisingly small. For decays to states with J > 3/2, these FBAs in 
general have more than one zero (for the states we have considered, there is one additional zero) . At leading order in 
HQET, the second zero does not exist. For decays to the ground state, it is known that the positions of the zeroes 
are modified in many scenarios that arise beyond the SM. This is also found to be true of the excited states, not 
surprisingly. Thus, it would be crucial to know the number of zeroes and their positions reliably (assuming that there 
will ever be sufficient statistics for the differential FBAs to be extracted with high precision). Our results indicate 
that the leading order predictions of HQET are misleading in this regard. 

The work presented here can be extended in several directions. We can use these form factors to study any number 
of polarization observables that can arise in these decays for leptons and baryons. Like the FBAs, these asymmetries 
are functions of the ratio of form factors and may therefore be less sensitive to the details of the form factor models. 
Such observables could provide more model-independent ways of extracting some of the Wilson coefficients. 

In comparing our rates with those of other authors, the main weakness of our model is highlighted, namely that 
the form factors predicted drop too rapidly as the momentum of the daughter baryon increases. This makes these 
form factors less reliable for nonleptonic decays such as Aj — >■ A(*V or the radiative decays Aj A(*'7, in which is 
small or zero. This is an inherent problem with models of this kind that use the harmonic oscillator basis. We need to 
emphasize that this is NOT an inherent flaw of all quark models, but of the particular choice of basis. Furthermore, 
this weakness is exacerbated by the larger range of q^ accessible in the rare decays of the Af,. It would be interesting 
to perform these calculations using wave functions computed from other bases, such as the Sturmian basis. The 
Sturmian basis leads to form factors whose kinematic dependence on appears in the form of multipoles. This has 
two significant advantages. First, it is common to expect form factors to have multipole dependence experimentally. 
Second, multipole forms would decrease significantly less rapidly as the momentum of the daughter baryon increases. 
These forms can therefore be expected to lead to larger branching fractions for many of the decay modes that we 
have considered. The semianalytic method developed here using the harmonic oscillator basis can be easily adapted 
for the Sturmian basis. 

The sensitivity of the rare decays we have studied to new physics that arises beyond the Standard Model can be 
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easily studied with the form factors that wc have calculated herein. The operator product structure easily incorporates 
scenarios from beyond the Standard Model, such as supersymmetric models, models with fourth generation quarks, 
and models with universal extra dimensions. These types of investigations can help provide constraints on any possible 
new physics scenarios. 



Acknowledgment 

Wc gratefully acknowledge the support of the Department of Physics, the College of Arts and Sciences, and the 
Office of Research at Florida State University. This research is supported by the U.S. Department of Energy under 
contracts DE-SC0002615. 



Appendix A: Semianalytic Treatment of Hadronic Matrix Elements 

The hadronic matrix elements (Eq. 90) can be written in the form 

h',h 



where the coefficients a;i(/i/) are the products of the normalization of the baryon state (the ^J2Eb^ of Eq. 77), the 

expansion coefficients (the 'qf'' of Eq. 79) and the various Clebsch-Gordan coefficients that appear in the parent 
(daughter) baryon wave function, and the indices hQi') contain all the relevant quantum numbers being summed over 

for the parent (daughter) baryon state. U^^l^"^"' is the spectator overlap given by 



and W^'.^J'^^^'^g^' is the interaction overlap, which takes the form 

Brriq p\ \ /■ . ..2 ,.2 .;,,4-i* / ~ 
'A (Aw 



x{sipA +p, s,0 I sTb I b{p, s^))C'^p yi^mM- (A3) 

Here, k\ = ax\'/ax'Q.x and aw = ^J{oi{ + a{,)/2. 

Prom the power series definition of the Laguerre polynomials, the integral over pp can be done analytically to give 

"p' rip 



r'=Or=0 P' P 



where Nni{a) are the normalization constants for the oscillator basis, c^^ are the coefficients in the power series 
definition of the Laguerre polynomials. 



X is the Gaussian integral, 



cfe =(-i)r rjn + k+l) 

""^ ^ ' r(n - r + l)r(A; + r + l)r(r + 1) ■ ^ ' 



Iirn;b)^l e-^^' = ^^^±^ , (A6) 
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Kp = appi /apiap, app> = ■sJ{oi^p + Q!^')/2, and the Kronecker deltas 5m^,mp are due to the orthogonality of the 

spherical harmonics. 

We can write the solid harmonics in Eq. A3 in terms of solid harmonics over k by making use of the addition 
theorem, 



A=0 fJ,\ = -\ 



where 



Iimil2m2 



Att{21 + 1)!! 



(2A + 1)!!(2I - 2A + i)!!'^VAi-Am-,.^3^i-Am-,.x W> 

(^ams I ^2m2 I ^limi) i 



and l^rn are the usual spherical harmonics. Eq. A3 therefore becomes 



exp 



Pa 

2m\ al w 



"a''a'So':AVa' , 



where 



•12 ] X't^x 



x{s{pA+P,Sq') I sTb I b(p,Sg))Cnx 



3^ama(A)- 



In its most general form, the quark current {s{p\ + p, Sqi) \ sTh \ b{p, Sq)) can be written 



(A7) 

(A8) 
(A9) 



(AlO) 



(All) 



{s{pK+P,Sq>) I sTh I b{p,Sq)) = E E ;^m/^' PA)3^£m, (^) • 

t=0 me=-e 



(A12) 



It was discussed in Section IV B 2 that the coefficients are complicated functions of pA • k. To see this, recall Eq. 
103, 



(s(PA +P, -1/2) I S7+756 I b{p, +1/2)) = — 

By use of Eq. A7, yim{p) can be written in terms of yim{k) as 

yiiip) = V3c{pA)+yoo{k)+yn{k), 

3^22 (P) = 



c'iPAY+yooik) + vwc{pA)+yii{k) + y22{k). 



Now {s{pa +p,—1/2) I S7+756 I +1/2)) becomes 

Mpa+P,-1/2)|s7+756|6(p,+1/2)) = -f{k,pA;x) 



STrc{c+l){pA)lyoo{k) 



+^/-{2c+i){pA)+yii{k) + \l^y22{k) 



167r, 
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where, 

f{k,pA;x) 



and 



= [(c + Ifpl + k^ + ml + 2(c + 1)paM'^'- 

In the above, we have used x = pA-k, where pa = pa/pa and k = k/k. Matching {s{pA+p, —1/2) | 57+ 756 | b{p, +1/2)) 
with the general form of the quark current (Eq. A12), we find that 

^-^hXsiOoC^'^a) = -VSTrc{c+l){pA)lf{k,pA;x), 

?-7+'75;1i(^'^a) = -W^(2c+ l)(pA)+/(fc,pA;a;), 



^m220,P\) = -\j^f{k,pA\x). 
The ^ coefficients can each be expanded in terms of Legendre polynomials as 



^tZS^^P^) = Y.C'Z,L{k,PA)PL{x), 



L=0 

where 

r+i 



CXVe^'PA) = (^^) cia;^X;(fc,pA)PL(a;) 



Using 

+L 



^^^"^^ = f im:) Yl^{pA)Yr.M{k) 
we can write the ^ coefficients in terms of Y^Mik), 



00 L 

eS(^'^^A) = E E C;CLM.(fc,??A)lLMjfc), (A13) 
L=OMi,=-I, 



where 



C''-lmeLM^{k,PA) = ( 2L + i ) ^LMip^KlmLikjA) 

= 2nYlj^^{pA) £\xer(e2Ak,PA)PL{x). (A14) 
The integral above is done numerically. The quark current can now be written as 

{sipA + P, s,. )\sTh\ b{p, s,)) = E E CrC^M. {k, Pa)Ylm, {k)yem, (k). (A15) 

line LMl 

We also note that the normalized Laguerre polynomials in Eq. All are functions of pa • k. These functions are 
expanded in terms of spherical harmonics as well, giving 



n L 



(^]=N^lW)J2 E ^nf{k,CPA)YLM{k), (A16) 
^ L=OM=-L 
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where 



+1 



dxL' 



i+h ( r_ 



Pl{x). 



Here, we have used Eq. 93. K^f^ can be determined analytically and the values for n < 1 are 

AO? = 47rFoo(PA), 



A 00 



47rro*o(PA) 



AIM ^ ^.l^P^Y^M- 



^' 3 ^2 

Substituting Eqs. A15 and A16 into Eq. All, we get 



L„iM„, J.-^ \ A LaM, 



XJ^AVv (k)YL^'M^, {k)YLM, (mm, (fc)i^L„M„ (fc)3^A^. (k). 



Recall that, 3^;^ = k^Yim- Making this substitution leads to 



imi LMl L^,M^, La Ma 
X{YX'IJ,^, I YL^,M^,YLMLYemeYLaMa I ^A/*^), 



(A'+€ + A + 2;4) 



where 



Sg,SgL^,Ma,LaMa 2n_ / ^jt-p " 4 t V^''' C tTa ^ 

-'^r;£m«I,Mi«;,,i;^mAUV"i'^Aj " '^^^ ^T;emtLMLy'^^P'^) 
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(A17) 



(A18) 



(A19) 
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(A21) 



is the momentum integral, which is done numerically. The bra-ket term is the angular integral over the spherical 
harmonics which can be handled analytically. We note here that the integral over the spherical harmonics provides a 
cut-off for the infinite sum over L in Eq. A13; this sum can be truncated &iL = \ + La+t + La' + A'. 



Appendix B: Analytic Form Factors 

Here we present the tensor form factors obtained from single-component quark model wave functions in the oscillator 
basis. The form factors for the vector and axial vector ciirrcnts have the same form as those published in [84]. These 
form factors were extracted using a nonrelativistic reduction of the quark current. The resulting form factors have the 
form of a polynomial in the daughter baryon momentum pA times a Gaussian. The polynomials in the form factors 
presented here are truncated at zeroth order in pA- 



Rare dileptonic decays of A(, in a quark model 



38 



1. 1/2"* 



Hi 
Hi 



Ih 



2^,2 



+ 1 



? 9 2 





= I 



H 



nioa 



9 9 9 



mbirisaly 



Ih = 



V "aa' 



3/2 



exp - 



pI 

2ml "aa' 



2. i/2r 



^4 



"aI M-M' 



-ff2 = Ih 



+ ■ 



'A I 7n ,x 



3ni 



H^ = /, 



Imlalal (al - al 
Gmbnisa^x' 



12mbmsa\y 



&mbmsa\y 

/3 / axax 



+ 



6mba\y 



3/2 



exp 



Pa 

2ml o^lx' 



3ni 



Gmga^y 



3. 1/2" 



Hi 

H2 



Ha 



Ih 
Ih 



H. = I 



mqttx I 3a^ — 2a^, 



H 



LH 



ax 
3mb 



ax 
3mb 



GmbmsO^ 
2m„ 



+ 



+ 



XX' 

ax 



ax 
2mb 



ax 2m. 
mgax (sal - 2a\, 
^mbmsO?xx' 



QmbmsOxx' 



2m^^ax 



2m„ 



ax 



ax 



mbUxx' ^TTT-b 



Rare dileptonic decays of A(, in a quark model 



39 



OL\a\' 



a 



5/2 



exp 



p\ 

2ml (^y 



4. 3/2" 



Ha 
He 



Ih 
Ih 
Ih 
Ih 
Ih 
Ih 



7mqa\a\, 3ax SrUq 



2 _ „ „.2 



mbmsafy 



2mbmsa\y 



mbnisaj^y 2mbmsaj^y 



■mbmsa\y 



a-x 
nib 



mbirisa^y 



^ 2ax 
nrib 

irn^ax 
rribaly 



rrisaxaly 
ax 



2mb 



iH 



1 / axax' 



5/2 



exp 



2ml «L' 



5. 3/2"* 



Hi 

H2 



Hi 
H, 

Ha 



Ih 
Ih 



2al 



2mbmsa\^. 
2mo 



5al 



Smq 
2mb 

5171, „ 



mb 



Vlmbms 2ms 



-TO' 



H^i = Ih 



5al 



2al, 



2mbm,sOL^ 



2 ' 5a2 _ 2a2 



rriq 
mb 



Ih 
Ih 
Ih 



6mq 



6to„ 



5al 



12TObTOs 



2mbmsaj^X' 



5al 



6ml 
mbaly 



+ 



2m„ 



6mbms 



2m„ 



5al 



mb 6mbms 



5ai 



mb 12mbms 



axax' 



a 



XX' 



7/2 



exp 



pI ' 
2r^A alx' . 



Rare dileptonic decays of A(, in a quark model 



40 



6. 5/2"* 



Hi 
He 



= I 



H 



= I 



3m„ Smf, 



+ 



•A' 



2 2 



H 



9 2 

K"a' 



nibmsOc 



mbrUsa^X' 



2m„ 



rub 



2mlal, 



mbrrisaly 



V2\ 



UA' 



7/2 



exp 



rub 



rrisalaly 



mbnisa^y mbnisai^y 



3m3 
rubaly 



rub 



p\ 

2ml "aa' 



Appendix C: Hadronic Tensors 

In this appendix we present the hadronic tensors we use in this work. Recall that the squared amplitude for 
dileptonic decays with unpolarized baryons is given by 



(CI) 



where L are the leptonic tensors. The most general Lorentz structure for each hadronic tensor H^'^ is 

H^-' = -a^g^''^ + Pi+Q'^Q" + ^i.Q^q" + pL+q^Q" + [iL^q^'q'' + h^s'"'"^Q„qp, (C2) 

where Q = pa^ +Pa is the total baryon 4-momentum and q = pa,, — Pa is the 4-momentum transfer. The coefficients 
QJ) /3±±) and 7 are written in terms of products of form factors. A generic coefficient A in the tensor Hj^'^ takes the 



form 
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Vl,4: 



%,1 



1 



40r2 



'/4,1 



40^5/2 
1 

'20r2 



(r2-2r(s + l) + (s-l)2)', 

(r2-2r(.§ + !) + (! -1)2)' 



1 



p (V^ - 1) - r5/2 + 2r3/2(s + 1) + 2r3 - r2(5s + 3) + 
4r(s - l)s - (s - 1)2 - (s - 1)'^^ . 



V2,l 
V2,3 
V2,4 



'/1, 2 
%,2 
'/4,2 



1 



40r2 



(^/F+l) (r2-2r(s + l) + (s-l)2)', 
^ (r-2^/F-s + l)''(r + 2^/F-s + l)^ 



80r5/2 
1 



^ ^6r3/2 + - 4r(s - 1) - 2Vr(s - 1) 
{s-lf^ (r-2^/F-s + l)^ 



''3,2 
'/3,4 



AA{^^ + ) 
'/1, 3 

'/2,3 



1 



40,5/2 (r2-2r(5 + l) + (^-l)2)% 

_l_(,_2VF-5 + l)'(r + 2^A^-s + l)^ 



(s- 1)2^ (r-2^/F-s + l)^ 
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rjff^^-^ = rjt^f^-+^ = -^{V^+l)^r'/'-2r'/\s + l)+r\s-l)-2r{rs'-l) + 

^AMP.-) = ^AAip.,) ^ ^2r3/^ + r' + " 1) " ^ 1)^) (r - 2^^ - s + l)\ 

<3 ^'^^-^ = Vtf'-^' = ^ (2r^/^ +r^+ 2V-r{s -!)-{§- if) (r - 2^? - J + l)^ 



= ^(r-2V^-s + ir(r + 2VF-a + l)\ 
j^y..) ^ _l_j'6r5/2_6r3/2s + 3r3 + r2(3-6s)+r(4s2-5s + l) + 



2v^(s-l)2-(s-l)3 . 



= ^(r-2VF-s + lf (r + 2Vf-s + l)^ 



2VF(s-l)2-(s-l)3^. 



Vir'--' = Vir--' = ^iV-r-l) {r^ - 2ris + 1) + - 1)^) 
vrr'--' = vir--'=-:^{V-r-l){r^-2ris+l)Hs-ir)\ 

VViP—) 



40r5/2 



' = ^.^^z^--) = -^(r'/' - r^'\28 + 3) + r^'^ + 3) - 3r3 + r\7s + 5) + 



r (-5s2 + 6s - 1) - v^(s - 1)2 + (s - 1)3 j , 
nX^s''--' = nV^''--' = (r - 2VF - a + 1)^ (r + 2Vf - ^ + l)^ 

vZ^'--^ = <2^^'--^ = - 45^:2 ( - + - 4r(s - 1) + 2V?(s - 1) + 

(s-l)2^ (r- + 2VF-s + l)', 

<r^^--^ = Vir--' = ( - 6-'^' + - 4r(s - 1) + 2Ms - 1) + 



{s-lf^ {r + 2^-s + lf . 
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Vtf'--^ = Vtf'--^ = 4^ (^/^ + 1) ir' - 2r(s + !) + {§- if)' , 

r (-5s2 + 6s - 1) + V^(s - 1)2 + (s - 1)3^ , 

rhf'--^ = Vtf'--^ = (e-'^' + 3^-' - - 1) - 2V-riS - 1) + 

(r--2VF-s + l)', 

<f ''--^ = Vtf'--^ = ^ (Sr^/^ + - 4r(s - 1) - 2^? - 1) + 
(s-l)2^ (r-2^/F-s + l)^ 



{r^-2r{s + l) + {s-lf) 



2 



^1,1 - ^1,1 - 20r2 

yA(7) _ AV{y) _ - 2r(g + 1) + (s - 1)^ 

^4,4 - 114,4 - 20r ' 

yA(7) _ Avi-r) _ - 2r(J + 1) + (1 ^ 1)^ 

Vi,4 - ri4,i - 

VA{^) _ AV{^) _ - 2r(g + 1) + (g - 1)' 

'?4,l - Vl,4 - 



[1] See, for example, the 2010 review by A. Ceccucci, Z. Ligeti and Y. Sakai in The Review of Particle Physics, K. Nakamura 
et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[2] See, for example, M. Neubert, Phys. Rcpt. 245, 259-396 (1994). 

[3] Sec, for example, N. Isgur, D. Scora, B. Grinstcin and M. B. Wise, Phys. Rev. D 39, 799 (1989); D. Scora and N. Isgur, 

Phys. Rev. D52, 2783-2812 (1995). 
[4] N. Barik and P. C. Dash, Phys. Rev. D 53, 1366 (1996). 
[5] See, for example, A. Radyushkin, Phys. Lett. B 271, 218 (1991). 

[6] See, for example, C. Aubin et al. [Fermilab Lattice and MILC and HPQCD Collaborations], Phys. Rev. Lett. 94, 011601 

(2005). 

[7] A. J. Buras and M. Munz, Phys. Rev. D 52, 186 (1995). 

[8] G. Buchalla, A. J. Buras and M. E. Lautcnbaclicr, Rev. Mod. Phys. 68, 1125 (1996). 

[9] A. F. Falk, M. Luke, M. J. Savage, Phys. Rev. D 49, 3367 (1994). 
[10] G. Greub, A. loannissian and D. Wyler, Phys. Lett. B 346, 149 (1995). 
[11] P. Gho, M. Misiak and D. Wyler, Phys. Rev. D 54, 3329 (1996). 
[12] W. Roberts and F. Ledroit, Phys. Rev. D 53, 3643 (1996). 
[13] W. Roberts, Phys. Rev. D 54, 863 (1996). 
[14] J. L. Hewctt and J. D. Wells, Phys. Rev. D 55, 5549 (1997). 
[15] T. M. Aliev, M. Savci and A. Ozpineci, Phys. Rev. D 56, 4260 (1997). 
[16] T. M. Aliev, G. S. Kim and Y. G. Kim, Phys. Rev. D 62, 14026 (2000). 
[17] T. M. Aliev, M. K. Gakmak, A. Ozpineci and M. Savci, Phys. Rev. D 64, 55007 (2001). 
[18] T. M. Aliev, A. Ozpineci and M. Savci, Phys. Lett. B 511, 49 (2001). 
[19] T. M. Aliev, M. K. Gakmak and M. Savci, Nuc. Phys. B 607, 305 (2001). 
[20] T. M. Aliev, A. Ozpineci, M. Savci and G. Yuce, Phys. Rev. D 66, 115006 (2002). 



Rare dileptonic decays of A(, in a quark model 



51 



D. Mclikhov, N. Nikitin and S. Simula, Phys. Lett. B 410, 290 (1997). 
D. Mclikhov, N. Nikitin and S. Simula, Phys. Lett. B 430, 332 (1998). 

D. Melikhov, N. Nikitin and S. Simula, Nuc. Phys. B (Proc. Suppl.) 75B, 97 (1999). 
T. Huang and Z.-H. Li, Phys. Rev. D 57, 1993 (1998). 

G. Burdman, Phys. Rev. D 57, 4254 (1998). 

T. Goto, Y. Okada and Y. Shimizu, Phys. Rev. D 58, 094006 (1998). 
C.-S. Huang, W. Liao and Q.-S. Yan, Phys. Rev. D 59, 011701 (1998). 
C.-S. Huang and Q.-S. Yan, Phys. Lett. B 442, 209 (1998). 

Q.-S. Yan, C.-S. Huang, W. Liao and S.-H. Zhu, Phys. Rev. D 62, 094023 (2000). 
C.-S. Huang, Nuc. Phys. B (Proc. Suppl.) 93, 73 (2001). 
C.-S. Huang, Nuc. Phys. B 657, 304 (2003). 
P. Colangelo et al, Eur. Phys. J. C 8, 81 (1999). 

P. Colangelo, F. De Fazio, R. Ferrandes and T. N. Phan, hep-ph/0709.2817v2 (2008). 
C.-K. Chua, X.-G. He and W.-S. Hou, Phys. Rev. D 60, 014003 (1999). 
A. Ah, P. Ball, L.T. Handoko and G. Hillcr, Phys. Rev. D 61, 074024 (2000). 
C.-H. V. Chang, D. Chang and W.-Y. Kcung, Phys. Rev. D 61, 53007 (2000). 

E. O. Itlan, Phys. Rev. D 61, 54001 (2000). 

M. Beneke, T. Feldmann, D. Seidel, Nuc. Phys. B 612, 25 (2001). 

F. Kruger and E. Lunghi, Phys. Rev. D 63, 14013 (2001). 
A. S. Safir, EPJdirect C15, 1 (2001). 

C.-H. Chen and C. Q. Ceng, Phys. Rev. D 63, 114025 (2001). 
C.-H. Chen and C. Q. Ceng, Phys. Rev. D 66, 14007 (2002). 
C.-H. Chen and C. Q. Ceng, Phys. Rev. D 66, 34006 (2002). 
C.-H. Chen and C. Q. Geng, Phys. Rev. D 66, 94018 (2002). 
M. Prank and S. Nie, Phys. Rev. D 66, 55001 (2002). 
A. Arda and M. Boz, Phys. Rev. D 66, 75012 (2002). 

S. R. Choudhury, N. Gaur, A. S. Cornell and G. C. Joshi, Phys. Rev. D 68, 54016 (2003). 
A. J. Buras, A. Poschenriedcr, M. Spranger, A. Weiler, hep-ph/0306158vl (2003). 

A. Ghinculov, T. Hurth, G. Isidori and Y.-P. Yao, Nuc. Phys. B 685, 351 (2004). 
V. Bashiry, hep-ph/0605061 (2006). 

F. De Fazio, Nuc. Phys. B (Proc. Suppl.) 174, 185 (2007). 

R. Ferrandes, hep-ph/0710.2212v2 (2007). 

M. J. Aslam, Y.-M. Wang and C.-D. Lu, Phys. Rev. D 78, 114032 (2008). 
Y. M. Wang, Y. Li and C.-D. Lu, hcp-ph/0804.0648 (2008). 
C. S. Huang, H. G. Yan, Phys. Rev. D 59, 114022 (1999). 
C.-H. Chen and C. Q. Geng , Phys. Rev. D 63, 114024 (2001). 
C.-H. Chen and C. Q. Geng , Phys. Rev. D 64, 074001 (2001). 
C.-H. Chen and C. Q. Geng , Phys. Lett. B 516, 327 (2001). 

H. Y. Cheng, C. Y. Cheung, G. L. Lin, Y. C. Lin, T. M. Yan, H. L. Yu, Phys. Rev. D 51, 1199 (1995). 

T. Mannel, S. Recksiegel, J. Phys. G 24, 979 (1998). 

X. G. He, T. Li, X. Q. Li, Y. M. Wang, Phys. Rev. D 74, 034026 (2006). 
F. Zolfagharpour and V. Bashiry, Nucl. Phys. B 796, 294 (2008). 
R. Ammar et a/.(CLEO Collaboration), Phys. Rev. Lett. 71, 674 (1993). 
M. Nakao et al. (Belle CoUab.), Phys. Rev. D 69, 112001 (2004). 

B. Aubert et al. (BaBar CoUab.), Phys. Rev. Lett. 103, 211802 (2009). 
B. Aubert et al. (BaBar CoUab.), Phys. Rev. D 70, 091105R (2004). 
H. Yang et al. (Belle CoUab.), Phys. Rev. Lett. 94, 111802 (2005). 
S. Nishida et al. (Belle CoUab.), Phys. Lett. B 610, 23 (2005). 
B. Aubert et al. (BaBar CoUab.), Phys. Rev. D 74, 031102R (2004). 
J. Wicht et al. (Belle CoUab.), Phys. Rev. Lett. 100, 121801 (2008). 
S. Chen et al. (CLEG CoUab.), Phys. Rev. Lett. 87, 251807 (2001). 
B. Aubert et al. (BaBar CoUab.), Phys. Rev. Lett. 97, 171803 (2006). 

A. Limosani et al. (Belle CoUab.), Phys. Rev. Lett. 103, 241801 (2009). 

E. Barberio et al. (Heavy Flavor Averaging Group), "Averages of b-hadron and c-hadron properties at the end of 2009," 
In preparation, and online update at http://www.slac.stanford.edu/xorg/hfag/. 
J.-T. Wei et al. (Belle CoUab.), Phys. Rev. Lett. 103, 171801 (2009). 

B. Aubert et al. (BaBar CoUab.), Phys. Rev. Lett. 102, 091803 (2009). 

T. Aaltonen et al.(CDF CoUab.), Phys. Rev. D 79, 011104 (2009); T. Aaltonen et aL(CDF CoUab.), hep-ex/1108.0695vl 
(2011). . 

M. Iwasaki et al. (Belle CoUab.), Phys. Rev. D 72, 092005 (2005). 
B. Aubert et al. (BaBar Collab.), Phys. Rev. Lett. 93, 081802 (2004). 

E. Barberio et al. (Heavy Flavor Averaging Group), "Averages of b-hadron and c-hadron properties at the end of 2007," 
arXiv:0808.1297 [hep-ex], and online update at http://www.slac.stanford.edu/xorg/hfag/. 
[81] A. Ishikawa et al. (Belle CoUab.), Phys. Rev. Lett. 96, 251801 (2006). 



Rare dileptonic decays of A(, in a quark model 



52 



[82] T. Aaltonen et al.(CDF CoUab.), hcp-ox/1107.3753vl (2011). 

[83] P. Koppenburg (LHCb Collaboration), LHC senuleptonic and radiative B decays program, The 7th International Conference 

on Hyperons, Charm and Beauty Hadrons (2006) 

[84] M. Pervin, W. Roberts and S. Capstick, Phys. Rev. C 72, 35201 (2005). 

[85] M. Pervin, W. Roberts and S. Capstick, Phys. Rev. C 74, 025205 (2006). 

[86] T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. B 355, 38 (1991). 

[87] W. Roberts and M. Pervin, Int. J. Mod. Phys. A, 23, 2817 (2008). 

[88] G. Crawford et aL(CLEO CoUaboration), Phys. Rev. Lett. 75, 624 (1995). 

[89] C. Amsler et al. (Particle Data Group), Physics Letters B667, 1 (2008). 

[90] See, for example, D. Jido, T. Sekihara, Y. Ikeda, T. Hyodo, Y. Kanada-En'yo, E. Oset, Nucl. Phys. A835, 59-66 (2010). 



